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THE ONSET OF INSTABILITY IN FIRST-ORDER SYSTEMS 


NICOLAS LERNER, TOAN NGUYEN, AND BENJAMIN TEXIER 

Abstract. We study the Cauchy problem for first-order quasi-Iinear systems of partial 
differential equations. When the spectrum of the initial principal symbol is not included in 
the real line, i.e., when hyperbolicity is violated at initial time, then the Cauchy problem is 
strongly unstable, in the sense of Hadamard. This phenomenon, which extends the linear 
Lax-Mizohata theorem, was explained by G. Metivier in [Remarks on the well-posedness 
of the nonlinear Cauchy problem, Contemp. Math. 2005]. In the present paper, we are 
interested in the transition from hyperbolicity to non-hyperbolicity, that is the limiting 
case where hyperbolicity holds at initial time, but is violated at positive times: under 
such an hypothesis, we generalize a recent work by N. Lerner, Y. Morimoto and C.-J. 
Xu [Instability of the Cauchy-Kovalevskaya solution for a class of non-linear systems, 
American J. Math. 2010] on complex scalar systems, as we prove that even a weak defect 
of hyperbolicity implies a strong Hadamard instability. Our examples include Burgers 
systems, Van der Waals gas dynamics, and Klein-Gordon-Zakharov systems. Our analysis 
relies on an approximation result for pseudo-differential flows, introduced by B. Texier in 
[Approximations of pseudo-differential flows, Indiana Univ. Math. J. 2016]. 
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1. Introduction 


We study well-posedness issues in Sobolev spaces for the Cauchy problem for first-order, 
quasi-linear systems of partial differential equations: 


( 1 . 1 ) 


E 

i<j<d 


dtu+ ^ Aj{t,x,u)dxjU = F{t,x,u), 


xWi X Miy into 


X X to the 


where t > 0, x € M“, u{t,x) G the maps Aj are smooth from 
space of X iV real matrices and F is smooth from 

We prove in this article a general ill-posedness result in Sobolev spaces for the Cauchy 


problem for (1.1), under an assumption of a weak defect of hyperbolicity that describes 


the transition from hyperbolicity to ellipticity. This extends recent results of G. Metivier 
|18| and N. Lerner, Y. Morimoto and C.-J. Xu [l2]. Here “well-posedness” is understood 
in the sense of Hadamard |3j, meaning existence and regularity of a flow; “hyperbolicity”, 
as discussed in Section EH means reality of the spectrum of the principal symbol, and 
“ellipticity” corresponds to existence of non-real eigenvalues for the principal symbol. 

We begin this introduction with a discussion of hyperbolicity and well-posedness (Sec¬ 


tion 1.1), then give three results: Theorem 1.2 describes ill-posedness of elliptic initial-value 


problems, while Theorems 1.3 and 1.6 are ill-posedness results for systems undergoing a 
transition from hyperbolicity to ellipticity. These results are illustrated in a series of exam¬ 
ples in Section 1.5 Our main assumption (Assumption |2.l| ) and main result (Theorem 2.2) 
are stated in Sections [Q and 12.21 

1.1. Hyperbolicity as a necessary condition for well-posedness. Lax-Mizohata theo¬ 
rems, named after Peter Lax and Sigeru Mizohata, state that well-posed non-characteristic 
initial-value problems for first-order systems are necessarily hyperbolic, meaning that all 
eigenvalues of the principal symbol are real. 

P. Lax’s original result [lOj is stated in a framework, for linear equations, i.e. such 
that Aj(t,x,u) = Aj{t,x). Lax uses a relatively strong definition of well-posedness that 
includes continuous dependence not only in the data, but also in a source. This allows 
him in particular to consider WKB approximate solutions; the proof of m shows that in 
the non-hyperbolic case, if the eigenvalues are separated, the norms of high-frequency 
WKB solutions grow faster than the norms of the datum and source, for any k. The 
separation assumption ensures that the eigenvalues are smooth, implying smoothness for 
the coefficients of the WKB cascade of equations. In the same (7°° framework for linear 
equations but without assuming spectral separation, S. Mizohata [16] proved that existence, 
uniqueness and continuous dependence on the data cannot hold in the non-hyperbolic case. 

Later S. Wakabayashi [26| and K. Yagdjian iniiEsi extended the analysis to the quasi- 
linear case, but it was only in 2005 that a precise description of the lack of regularity of 
the flow was given, by Metivier: Theorem 3.2 in [T8| states that in the case that the Aj 
are analytic, under the assumption that for some fixed vector G and some frequency 
^0 g -^d principal symbol Yli<j<d^ji'^^)^j hyperbolic, some analytical data 


uniquely generate analytical solutions, but the corresponding flow for (1.1) is not Holder 
continuous from high Sobolev norms to L^, locally around a Canchy-Kovalevskaya solution 
issued from the constant datum 
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Metivier’s result is a long-time Cauchy-Kovalevskaya result. Without loss of generality, 
assume indeed that vP = 0. Then Theorem 3.2 in [18] states that data that are small in high 
norms may generate solutions that are instantaneously large in low norms. In this view, 
assume in (1.1) the hyperbolic ansatz: u{t,x) = sv{t/e,x/s), where e > 0. Setting F = 0 
for simplicity, and t = tje, y = xfe, the equation in v is 


( 1 . 2 ) 


drV + ^ ^j(eT, ey, ev)dy^v = 0. 
l<3<d 


If all fluxes Aj are analytic in their arguments, the Cauchy-Kovalevskaya theorem ensures 
the existence and uniqueness of a solution v issued from an analytic datum v{0,x), over a 
time interval 0(1) in the fast variable r. What’s more, by regularity of the coefficients of 
(1.2) with respect to e, the solution v stays close, in analytical semi-norms, to the solution 


w of the constant-coefficient system 

drW + ^ Aj{0)dy.w = 0 

over time intervals 0(1). By Assumption on Aj(0), the Fourier transform of w in 

the spectral direction grows like for some O > 0. This implies a similar growth 

for v{t/e, ^®), and in turn a growth in for u{t, but only on time intervals 0(e), 

due to the initial rescaling in time. The content of Metivier’s result is therefore to show that 


the solution v to (1.2) exists, and the growth persists, over “long” time intervals 0(| loge|). 


so that the exponential amplification is effective. 

In the scalar complex case, the results of N. Lerner, Y. Morimoto and C.-J. Xu |12) 
extended the analysis of Metivier to the situation where the symbol is initially hyperbolic, 
but hyperbolicity is instantaneously lost, in the sense that a characteristic root is real at 
t = 0, but leaves the real line at positive times. The main result of m states that such a 
weak defect of hyperbolicity implies a strong form of ill-posedness; the analysis is based on 
representations of solutions by the method of characteristics, following mi- This argument 
does not carry over to systems, even in the case of a diagonal principal symbol, if the 
components of the solution are coupled through the lower-order term F{u). 


Our goal in this article is to extend the instability results of [T2| on complex scalar 
equations to the case of quasi-linear first-order systems (1.1). In the process, we recover a 
version of the results of 


, with a method of proof that does not rely on analyticity. 

1.2. On the local character of our assumptions and results. Our assumptions are 
local in nature. They bear on the germ, at a given point {to,xo,^Q) G M+ x M'’* x MP, 
representing time, position, and frequency, of the principal symbol evaluated at a given 
reference solution. Under these local assumptions, we prove local instabilities, which extend 
the aforementioned Lax-Mizohata theorems, and which roughly say that there are no local 
solutions possessing a minimal smoothness with initial data taking values locally in an 
elliptic region. These local instabilities are independent of the global properties of the 


instance the compressible Euler equations (1.16) introduced in Section 1.5 


system (1.1). In particular, the system (1.1) may have formal conserved quantities; see for 
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1.3. Transition from hyperbolicity to ellipticity. Our starting point is to assume that 
there exists a local smooth solution cj) to 0 with a large Sobolev regularity: 

(1.3) <PGC°°i[0,n],H^^{U)), 

for some Tq > 0, some open set U C and some Sobolev regularity index si = l+d/2-\-S2, 
where S 2 > 0 is large enough, depending on the parameters in our probleirQ If the matrices 
Aj and the source F depend analytically on {t, x, u), then we can choose (/> to be a Cauchy- 
Kovalevskaya solution. However, we do not use analyticity in the rest of the paper. 

The linearized principal symbol at 4> is 

(1.4) A{t,x,^):= ^ ^jAj{t,x,(l){t,x)). 


The upcoming ill-posedness results are based on readily verifiable conditions bearing on the 
jet at f = 0 of the characteristic polynomial P of the principal symbol: 


(1.5) 


P{t, X, A) := det (A Id —A{t, x, ^)). 


Most of these conditions are stable under perturbations of the principal symbol, and all can 
be expressed in terms of the fluxes Aj and the initial datum ^(0). In particular, it is of key 
importance that the verification of these conditions does not require any knowledge of the 
behavior of the reference solution (p at positive times. 

Also, it should be mentioned that our hypotheses do not require the computation of 


eigenvalues and are expressed explicitly in terms of derivatives of P given by (1.5) at initial 
time. 


1.3.1. Hadamard instability. If (1.1) does possess a flow, how regular can we reasonably 


expect it to be? A good reference point is the regularity of the flow generated by a symmetric 
system. If for all j and all u, the matrices Aj(u) are symmetric, then local-in-time solutions 
to the initial-value problem for (1.1) exist and are unique in H®, for s > 1 + d/2 [3l[71[9]; 


moreover, given a ball B}{‘< (0,4?) C H^, there is an associated existence time T > 0. The flow 
is Lipschitz Bh‘>{0,R) Ci —)• L°°{[0,T], H^), continuous Bhs(0,R) —)• L°°{[0,T], H^), 

but not uniformly continuous Bh»{0,R) —)> L°°([0, T], H*) in general [7]. Micro-locally 
symmetrizable systems also enjoy these properties |19j . 


Accordingly, ill-posedness will be understood as follows: 


Definition 1.1. We will say that the initial-value problem for the system (1.1) is ill-posed 


in the vicinity of the reference solution f satisfying (1.3), if for some xq G U, given any 
parameters m,a,6 > 0, T such that 


( 1 . 6 ) 


m G 


^ < a < 1, 


B{xq,5) cU, 0 < T < To, 


where U and Tq are as in (1.3), there is no neighborhood U of (j){0) in such that, 

for all u{0) G Z?, the system {1.1) has a solution u £ L°°{[0,T],W^’°^{B{xo,6))) issued from 


^We use regularity of p in p articu lar in the construction of the local solution operator; see Appendix [d| 
specifically the proof of Lemma D.2 in which go is the order of a Taylor expansion involving (f>. 
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u{0) which satisfies 
(1.7) 


sup 

UQ^Vl 

0<t<T 


\u{t) - 4>it)\\w^’<=°{B{xo,S)) 


Iko - <('(0)11 


< + 00 . 




Thus (1.1) is ill-posed near the reference solution if either data arbitrarily close to 0(0) 


fail to generate trajectories, corresponding to absence of a solution, or if trajectories issued 
close to 0(0) deviate from 0, corresponding to absence of Holder continuity for the solution 
operator. In the latter case, we note that: 

• the deviation is relative to the initial closeness, so that 0 is unstable in the sense of 
Hadamard, not in the sense of Lyapunov; 

• the deviation is instantaneous: T is arbitrarily small; it is localized: 6 is arbitrarily 
small. 

• The initial closeness is measured in a strong norm, where m is arbitrarily larg^ 

while the deviation is measured in a weaker norm, defined as |f|i^/i,oo = 

|/|l- + |V,/|loo. 

In our proofs of ill-posedness in the sense of Definition |l.l[ we will always assume existence 
of a solution issued from a small perturbation of 0(0), and proceed to disprove (1.7). 


Note that the flows of ill-posed problems in the sense of Definition o exhibit a lack of 
Holder continuity. F. John introduced in [6] a notion of “well-behaved” problem, weaker 
than well-posedness. In well-behaved problems, Cauchy data generate unique solutions, 
and, in restriction to balls in the W^'°° topology, for some integer M, the flow is Holder 
continuous in appropriate norms. The notions introduced in [6] were developed in the article 
[T] by H. Bahouri, who used sharp Carleman estimates. 

The restriction to a > 1/2 in Definition is technical. Precisely, it comes from the 
fact that we prove ill-posedness by disproving as indicated above. This gives weak 

bounds on the solution, which we use to bound the nonlinear terms. Consider nonlinear 
terms in (1.1) which are controlled by I'o-homogeneous terms in u, with > 2, that is 
such that duAj = and duF = These bounds hold if, for instance, 

Aj{u)dx.u = u^^~^dxU and F{u) = using scalar notation. Then, the proof of our 


general result (Theorem 2.2) shows ill-posedness with a > l/io- (See indeed Lemma 3.16 
and its proof, and note the constraint 2K' > K which appears at the end of the proof in 


Section 3.15 


Finally, we point out that Definition 1.1 describes only the behavior of solutions which 
belong to This excludes in particular shocks, which are expected to form in finite 


time for systems (1.1), even in the case of smooth data. Shocks with jump across elliptic 


zones could exhibit some stability properties. 

1.3.2. Initial ellipticity. Our hrst result states that the ellipticity condition 
(1.8) P(0,wo) = 0, Wo = (a:o,?o,Ao) G H X (M'^Z {0}) X (C\M), 


where P is the characteristic polynomial defined in (1.5), implies ill-posedness: 


^That is, the only restriction on m is the Sobolev regularity of 0 : we need, in particular, m < si for (1.71 


to make sense. 
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Figure 1. In Theorem |1.2[ corresponding to £ = 0 in Assumption 2.1 the 


principal symbol at (0,xo,^o) has non-real eigenvalues Ao,Ao. These may 
correspond to coalescing points in the spectrum, for near (0,xo,Co)- 


Theorem 1.2. Under the ellipticity condition (1.8), the Cauchy problem for system 
is ill-posed in the vicinity of the reference solution in the sense of Definition [Tt} 

Theorem |1.2| (proved in Section]^ states that hyperbolicity is a necessary condition for the 
well-posedness of the initial-value problem (1.1), and partially recovers Metivier’s resul10 
An analogue to Theorem |1.2| in the high-frequency regime is given in based on [ 23 ] just 


like our proof of Theorem 1.2; the main result of [T3| precisely describes how resonances may 


induce local defects of hyperbolicity in strongly perturbed semi-linear hyperbolic systems, 
and thus destabilize WKB solutions. 


1.3.3. Non semi-simple defect of hyperbolicity. We now turn to situations in which the 
initial principal symbol is hyperbolic: 

(1.9) P(0, x,^,A) = 0 implies A G M, for all (x, G 17 x (M'^ \ {0}), 
and aim to describe situations in which some roots of P are non-real for t > 0. Let 

r := {w = (x,^,A) G 1/ X (M‘^\{0}) X M, P{0,uj) =0}, 

By reality of the coefficients of P, non-real roots occur in conjugate pairs. In particular, 
eigenvalues must coalesce at t = 0 if we are to observe non-real eigenvalues for t > 0. 

Let then cjq G F, such that 

( 1 . 10 ) dxP{0,uo) = 0 , dlP{0,iOo) / 0 . 

The eigenvalue Aq of A(0, xo,f,o) thus has multiplicity exactly two. Assume in addition that 
loq satisfies condition 

(1.11) {dlPdtP){o,coo)>o. 

^Theorem 3.2 in [TH] shows not only instability, but also existence and uniqueness, under assumption of 
analyticity for the fluxes, the source and the initial data. 
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Figure 2. In Theorem |1.3[ corresponding to ^ = 1/2 in Assumption 2.1 


bifurcation occurs at (0, xq, Co) in the spectrum of the principal symbol. The 
eigenvalues are not time differentiable. The arrows indicate the direction of 
time. 


The eigenvalues are continuouQ implying that condition ( |1.11[ ) is open, meaning that if 
it holds at loq, then it holds at any nearby lj in T. 


Theorem 1.3. Assume that conditions (l.lO)-(l.ll) hold for some wq G T, and that the 
other eigenvalues of A{0,xo, Co) ore simple. Then the Cauchy problem for system 0 is 
ill-posed in the vicinity of the reference solution 0, in the sense of Definition 0 


The conditions (1.10)-(1.11) are relevant, and, as far as we know, new, also in the linear 


case. 


Van der Waals systems and Klein-Gordon-Zakharov systems illustrate Theorem 1.3 


see 


Sections 1.5, 7.3 and 7.4 


The proof of Theorem 1.3, given in Section]^ reveals that under ( 1.10[ )-( [1.11[ ), the eigen¬ 
values that coalesce at t = 0 branch out of the real axis. The branching time is typically not 
identically equal to t = 0 around {xo,Co)] for {x,C) close to (xq, Co)A^ is equal to t^(x,C) > 0, 
with a smooth transition function t*. At (C(x, C), x, C) the branching eigenvalues are not 
time-differentiable, in particular not semi-simple. Details are given in Section 5.1 in the 
proof of Theorem |1.3[ Figure 3 pictures the typical shape of the transition function. The 
elliptic domain is {t > t*}, and the hyperbolic domain is {t < t*}. 


Under the assumptions of Theorem |1.3| and assuming analyticity of the coefficients, 
B. Morisse m proves existence and uniqueness in addition to instability in Gevrey spaces, 
further extending G. Metivier’s analysis |18j . 


1.3.4. Semi-simple defect of hyperbolicity. Time-differentiable defects of hyperbolicity of 
size two can be simply characterized in terms of P : 


Proposition 1.4. Let P{t,x,C, X) be the characteristic polynomial (1.5) of the principal 
symbol A{t,x,C) (1-4). We assume initial hyperbolicity (1.9). Let oj = (x,^,A) G F. If 


^By continuity of A and Rouche’s theorem; see [8] or |25| . 
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Figure 3. In Theorem 1.3, the transition occurs at t = > 0, near (xo,^o)- 


d\P{0,uj) = 0 ^ (9|P(0,a;), then for the branches A of eigenvalues of A which coalesce at 
(0, X, f), there holds 

is differentiable ot t = Ol f dtP{0,uj) = 0, 

a \ / n r { /C 52 D^n . ,^^2 ^ do2 ; 


QmdtX{0,x,f) 7 ^ Oj {dfxP{0,^^)) < idfPdxP){0,uj). 

Proof. We assume X £ C^. The proof in the general case is postponed to Appendix]^ For 
t in a neighborhood of 0, there holds P(t, x, A(t, x, ^)) = 0. Differentiating with respect 
to t, we find 

dtP{0, u) + 9iA(0, X, f)dxP{Q, uj) = 0. 

Since A(0,x,^) is real-valued, by reality of P, the derivatives dtP and d\P are real. If we 
assume (9tA(0, x, / 0, then there holds dtP{oj) = d\P{oj) = 0. Differentiating again 
with respect to t, we find 

(1.13) 52p(0,u;) + 29iA(0,x,09u^’(0,w) + {dtXfd,x,i)fdlPfd,uj) = 0. 


Equation (1.13), a second-order polynomial equation in 9tA(0,x,^), has non-real roots if 
and only if the second condition in the right-hand side of (1.12) holds. □ 


We now examine the situation in which a double and semi-simple eigenvalue Aq belongs 
to a branch A of double and semi-simple eigenvalues at t = 0, which all satisfy conditions 


(1.12), that is: 


Hypothesis 1.5. For some ojq = (xo,^0)-^o) £ T satisfying (1.10) and (1.12), and such 


that Ao is a semi-simple eigenvalue of A{0,xo,fo), for all u = (x,^. A) in a neighborhood of 
Wo m F, there holds 

dxP{Q,uj) = 9tP(0,w) = 0, 
and X is a semi-simple eigenvalue o/A(0,x,^). 

Semi-simplicity of an eigenvalue means simpleness as a root of the minimal polynomial. 

Condition {df^P{uj))‘^ < {d]‘Pd\P){Ljj) is open; in particular if it holds at wq G F, it holds 
at all nearby w G F. Thus under Hypothesis 1.5 conditions (1.10) and (1.12) hold in a 
neighborhood of wq in T. 
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Figure 4. In Theorem |1.6[ corresponding to I" = 1 in Assumption |2.1[ a 
bifurcation occurs at (0, x, in the spectrum of the principal symbol, for 
all (x,^) near (xo,Co)- The eigenvalues are time-differentiable. The arrows 
indicate the direction of time. 




t 


f = 1 


elliptic 


(a^oTo) 

hyperbolic 

{x, 


Figure 5. 


In Theorem 1.6, the transition occurs at t 


0, uniformly near (xo,Co)- 


Theorem 1.6. Assume that Hypothesis \l.f)\ holds, and that the other eigenvalues o/A(0, xq, ^o) 
are simple. Then the Cauchy problem for system a is ill-posed in the vicinity of the ref¬ 
erence solution (j), in the sense of Definition a 


An analogue to Theorem 1.6 in the high-frequency regime is the result of Y. Lu |T3], in 
which it is shown how higher-order resonances, not present in the data, may destabilize 
precise WKB solutions. 


Theorem 1.6 is illustrated by the Burgers systems of Sections |1.5| and 7.1 


1.4. Remarks. Taken together, our results assert that, for principal symbols with eigen¬ 
values of multiplicity at most two, if one of 


(a) condition (1.8), 

(b) conditions ( 1.10[ )-( [m| ), 

(c) Hypothesis 


1.5 


holds, then ill-posedness ensues. 
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We note that condition (1.11) is stable by perturbation, and that conditions (1.10)-(1.11) 


are generically necessary and sufficient for occurence of non-real eigenvalues in symbols that 
are initially hyperbolic. Indeed: 

• non-real eigenvalues may occur only if the initial principal symbol has double eigen¬ 


values, implying necessity of conditi on ( 1.10), and 
as shown by the proof of Theorem 


condition (1.11) implies real eigenva 


1.3 


the opposite sign {d\PdtP){t),ijJo) < 0 i 


m 


ues for small t > Q. 


Here generically means that the above discussion leaves out the degenerate case dtP{0, wq) = 

0 . 


We consider the case dtP{0,oJo) = 0 in Theorem 1.6 
significant gap between (b) and (c), the assumptions of Theorems 1.3 and 1.6 


Note however that there is a 

Indeed, 


while condition dfP = 0 in Hypothesis 1.5 lies at the boundary of the case considered in 


Theorem 1.3 Hypothesis 1.5 describes a situation which is quite degenerate, since we ask 
for the closed conditions d\P = 0, dtP = 0 (and also for semi-simplicity) to hold on a whole 
branch of eigenvalues near Aq- 

Non-semi-simple eigenvalues are typically not differentiable at the coalescing point, the 
canonical example being 


(1.14) 


0 


with a = 1. The proof of Theorem 1.3 shows that the principal symbol at (t,a: 0 )^o) can be 
reduced to (1.14), with a = 1 and a negative sign, implying non-real, and non-differentiable 
eigenvalues. 

By constrast, semi-simple eigenvalues admit one-sided directional derivatives (see for 
instance Chapter 2 of T. Kato’s treatise [S], or j22L I25jb In particular, there is some 


redundancy in our assumptions of semi-simplicity and condition (1.12). 


We finally observe that our analysis extends somewhat beyond the framework of The- 
and |1.6[ Consider for instance in one space dimension a smooth principal 


orems 
symbol of form 


0 


x^t — + Pa{x) 0 


Q ^ , with eigenvalues X± = ±^(x^t — t‘^ + t^a{x)) 


1/2 


with a G M, so that the eigenvalues are time-differentiable only at x = 0 : conditions 
hold only at x = 0. Semi-simplicity does not hold at (t, x) = (0,0). Condition ( |1.11[ ) does 
not hold at {t, x, A) = (0, 0, 0). However, by the implicit function theorem, eigenvalues cross 
at (s(x), x) for a smooth s with s{x) = x'^ + O(x^). By inspection, condition ( |1.11[ ) holds at 
(s(x),x). Since x is arbitrarily small. Theorem 1.3 applies, yielding instability. 


1.5. Examples. Burgers systems. Our first example is the family of Burgers-type systems 
in one space dimension 

j- -6(u)^U2\ o fuf 

.U 2 ) \U2 ui / ^ 


(1.15) 


dt 


= F{u), 
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in which b > 0 and F = (Fi, F 2 ) G are smooth. When h{u) is not constant, the instability 
result for scalar equations in m does not directly apply. We show in particular that if F 2 is 
not initially zero, then Theorem 1 1. 6| yields ill-posedness for the Cauchy problem for (1.15). 


Under the same condition F 2 {t = 0) ^ 0, Theorem 1.6 also applies to two-dimensional 
systems 

uidx^ -h{ufu2{dx^ + 8x2) 

1^2(^xi T 8x2) Ui8xi 


8 tu + 


u = F{u). 


Details are given in Section 7.1 and|7.2 


(1.16) 


Van der Waals gas dynamics. Our results also apply to the one-dimensional, isentropic 
Euler equations in Lagrangian coordinates 

8tui + 8 xU2 = 0 , 

8tU2 + 8xp{ui) = 0 , 

with a Van der Waals equations of state, for which there holds p'{ui) < 0, for some ui G M. 
We prove that if {cfii, (j) 2 ) is a smooth solution such that, for some xq G M, there holds 

(i)p'((f>i(xo, 0 )) < 0, or (u)p'((/)i( 0 ,xo)} = 0, p"((/)i( 0 , xo)) 8 x(/) 2 ( 0 , xq) > 0, 

then the initial-value problem for ( 1.16[ ) is ill-posed in any neighborhood of cj). Condition (i) 
is an ellipticity assumption (under which Theorem |1.2| applies), and condition (ii) is an open 
condition on the boundary of the domain of hyperbolicity (under which Theorem |1.3| ap¬ 
plies). System (1.16) possesses the formal conserved quantity / {\v{t,x)\‘^ + 2P{u{t,x))) dx, 


where P' = p. As briefly discussed in Section 1.2, the instabilities we put in evidence are 
local and do not preclude nor are contradicted by global stability properties of the system, 


such as formal conservation laws. This example is developed in Section 7.3 

Klein-Gordon-Zakharov systems. Our last class of examples is given by the following 
one-dimensional Klein-Gordon systems coupled to wave equations with Zakharov-type non- 
linearities: 

a 0 
0 0 


(1.17) 


8 t 


+ 8 x 


+ 


8 x 


n 

m 


= (n-Fl) 


8 t 


n 

m 


+ c 8 x 


m 

n 




8 x 


= 8 x 


u 


V 

—u 

0 

-I- 


The linear differential operator in (n, v) in the subsystem in (u, v) is a Klein-Gordon op¬ 
erator, with critical frequency scaled to 1. The linear differential operator in {n,m) in the 
subsystem in (re, rre) is a wave operator, with acoustic velocity c. The source in 8 x{u^ + v'^) 
is similar to the nonlinearity in the Zakharov equation |20l I23j . Systems of the form (1.17), 
with a = 0, are used to describe laser-matter interactions; in the high-frequency limit, 
they can be formally derived from the Maxwell-Euler equations [2]. We consider the case 
|c| < 1, corresponding to the physical situation of an acoustic velocity being smaller than 
the characteristic Klein-Gordon frequencies. 

It was shown in [2] that, for a = 0, system (1.17) is conjugated via a non-linear change 
of variables to a semi-linear system, implying in particular well-posedness in i7^(M), for 
s > 1/2. 
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Here we show that if a 7 ^ 0 and cj) = (it, v, n, m) is a smooth solution such that 


u( 0 ,xo) = 0 , i;( 0 ,xo) =-, acdxuiQjXo) > 0 , for some xg G 

2 a 


then Theorem 1.3 applies and the Cauchy problem for (1.17) is ill-posed in the vicinity of (j). 
This situations is analogous to the Turing instability, in which 0 is a stable equilibrium for 
both ordinary differential equations X' = AX and X' = BX, but not for X' = {A + B)X. 


We come back to this example in detail in Section 7.4 


2. Main assumption and result 


Theorems 1.2, 1.3 and 1.6 can all be cast in the same framework, which we now present. 


2.1. Bounds for the symbolic flowf of the principal symbol. 


2 . 1 . 1 . Degeneracy index and associated parameters. Let i £ {0, 1 / 2 , 1 }q Associated with i, 
dehne 

1 , iff = 0 , 

2/3, if £ = 1 / 2 , and C = 


h = 


1 


1 + ■ 


= 


1 / 2 , if£ = l. 

Parameters h and define our time, space and frequency scales 


0, if£E{0,l}, 

1/3, if £ = 1/2. 


2.1.2. The time transition function and the elliptic domain. Introduce a time transition 
function t* such that 

• if£ = 0 or£=l, then £* = 0 , 

• if £ = 1 / 2 , then t* depends smoothly on (x,^) and singularly on e, and is slowly varying 
in X, in the sense that there holds for some smooth function 0 * : 

( 2 . 1 ) 

t*(e,x ,0 = e with > 0 , 6 '*( 0 ,Co) = 0 , ^ 0 ) = 0 , if £ = 1 / 2 . 

Define then the elliptic domaiiij^by 

(2.2) T> := {(r;i,x,0, U{x,f)<T<t<T{e), |x| < 5, 

for some ^0 £ \ { 0 }, some 5 > 0 , with 

(2.3) r(e)^+i = r*|loge|, T* > 0. 


1.6 


^Theo rem 1.2 corresponds to the case 1 = 0, while Theorem 1.3 corresponds to the case £ = and 
Theorem 


to £ — 1. 


®The elliptic domains corresponding to Theorems 1.3 and 1.6 are pictured on Figures 3 and 5. 
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2.1.3. The rescaled and advected principal symbol. We consider a reference solution (j) satis¬ 
fying (1.3), and the associated principal symbol (1.4). The rescaled and advected principal 
symbol i^ 

(2.4) A(e,f,x,0 := (Q{A - {eH, xq + x, 0 ), 


(2.5) 


dt 


,1-hf, .. I (t,xo + £^ 


(O,x,0 = 


for some xq G IR“, with 0 < h < 1 as in ( |2.3[ ), where, in domain P : 

• the symbol // = /i(f,x,^) is real and smooth, and Q{t,x,^) G £,NxN jg g];nooth and 
pointwise invertible, 

• the bicharacteristics (x*,^*) solve 

-d^li 

We assume that the symbol is block diagonal, and for the blocks A.^j of A* consider 
either the bound 

(2.6) £^~^\d'^d^Ai,j\ < Cap < oo, for some Cap > 0, in V, uniformly in e, 

or the block structure 

(2.7) 




a:* 


-h-i 


Ai,j — 


0 

£^~^Ai,j2l 


.h-1 


^*112 

0 


with \d^d^A *ji 2 \ + I^x^f^*i 2 i| < Cap for some Cap > 0, in V, uniformly in e. 

If ^ = 0, then h = 1. As a consequence, in the block diagonalization of A* all blocks 
satisfy (| 2 . 6 [), by the assumed smoothness of the components of A*. 


If i = 112, we assume that some block of A* satishes (2.7), and the other blocks of A* 
satisfy (| 2 . 6 |). 


If £ = 1, we assume that all blocks of A* satisfy (2.6). 

2.1.4. Symbolic flow and growth functions. The symbolic flow S = 5(r;e,t,x,^) of A* is 
dehned as the solution to the family of linear ordinary differential equations 

( 2 . 8 ) dtS + i£^~^Ai,{£,t,x,f,)S = 0, 5(r;r)=Id. 

In the above Section 2.1.3[ we assumed that A* is block diagonal. Accordingly, the solution 
S to (2.8) is block diagonal, with blocks S'(i), 5 ( 2 )) ■ ■ • 

Let 7 ^(x,^) be two continuous functions defined on {|x| < <5, — ^o| < snch that 

7"(0,^o) = 7’^(0,lo), and 


(2.9) 


(r; t, X, f) := exp (^ 7 =^ {x,f){{t- tflx, 0)++^ - (t - tflx, 0)+“^^)) , 


where t^ := max(f, 0) and the time transition function f* is defined in ([ 23 . We understand 
e..^± as growth functionrX measuring how fast the solution 5 to (2.8) is growing, as seen in 


Assumption 2.1 below. The associated 7 * are rates of growth. 


'In the elliptic case, corresponding to Theorem 


1.2 


we have = 0, /i = 1, Q = Id, /r = 0, so that 
= (x,(), and then A* is simply a:,^) = A(et,x,^). 

®In the elliptic case, we have £ = 0, = 0, so that the growth functions are simply e ± = 
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2.1.5. Bounds. We postulate bounds for S in the elliptic domain T> in terms of the growth 
functions e.^± : 

Assumption 2.1. For some (xo,^o) ^ U x \ {0}), some £, and as above, 

any T* > 0, for some (5 > 0, for symbol A* satisfying the structural assumptions of Section 


2.1.3, for some eo > 0, there holds for all 0 < e < eq, for the solution S to {2.8): 


• the lower bound, for some smooth family of unitary vectors e{x) G C^, for |x| < <5 : 

(2.10) e"^e.^-(0;r(e),x,^o) < | ^(O; r(e), x, Co)e(x) | , 

• the upper bound for the jth diagonal block Sj of S, for (r, t, x, G V : 

1 


( 2 . 11 ) 


|5'(j)(r;t,x,0l 


1 


{r]t,x,^). 


In (2.10), notation a < b, where a and b are functions of (e,r,t,x,^), is used to mean 


existence of a uniform bound 

(2.12) a(e, r, t, x, ^) < C\ lne|^ 6(e, r, t, x,^), 

where C > 0 and C" > 0 do not depend on {e,T,t,x,f). This means in particular that 
powers of | lne| play the role of constants in our analysis. They are indeed destined to be 
absorbed by arbitrarily small powers of e. 


In (2.11), we use notation < for matrices. Here we mean block-wise inequalities “modulo 


constants”, in the sense of (2.12). That is, in (2.11) we assume that the j-th diagonal block 
of S itself has a block structure, with the top left block being bounded entry-wise by e..^+, 
the top right block being bounded entry-wise by e“^e..y+, etc. 


We further comment on Assumption 2.1 in Section 2.3 


2.2. Hadamard instability. The non-linear information contained in Assumption 2.1 


on 


quasi-linear system (1.1): 


the symbolic flow of the principal symbol (1.4) transposes into an instability result for the 


Theorem 2.2. Under Assumption 2.1, system (O) is ill-posed in the vicinity of the refer¬ 
ence solution (f), in the sense of Definition [Tt} 


Theorem 2.2 states that either there exists no solution map, or the solution map fails to 


enjoy any Holder-type continuity estimates. The proof of Theorem |2.2| is given in Section 
Key ideas in the proof are sketched in Section |2.4[ 


Theorems 1.2, 1.3 and 11.61 all follow from Theorem 2.2 


2.3. Comments on Assumption 2.1, Our main assumption is flexible enough to cover the 


three different situations described in Theorems 1.2, 1.3 and 1.6 Before further commenting 


on its ingredients in Section 2.3.1 and its verification in Section 2.3.2 we point out two key 
features: 


• Assumption 2.1 is nonlinear. It bears on the whole system (1.1), not just the principal 
symbol. For instance, instability occurs for the Burgers systems of Section 1.5 under a 
condition bearing on the nonlinear term F. 

• Assumption |2 .1 1 is finite-dimensional, in the sense that it postulates bounds for solutions 
to ordinary differential equations in a finite-dimensional setting. These are turned into 
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bounds for the solutions to the partial differential equations via Theorem D.3 An informal 
discussion of the role of Theorem D.3|is given in Section 2.4 


2.3.1. On the ingredients of Assumption [2T| 

• Our localization constraints in (x, .^) G respect the uncertainty principle. Indeed, 

we localize spatially in a box of size ~ We localize in frequency in a box of size 

~ around but then in the proof use highly-oscillating data and an e^-semi-classical 
quantization, so that it is which belongs in a box of size around meaning a 
frequency localization in a box of size If ^ = 0 or .£ = 1, then = 0. The area of 

the (x,^)-box is then £ 2 d(i- 2 /i) ^ since /i = 1 or h = 1/2. If £ = 1/2, then h = 2/3 and 
( = 1/3. The area of the (x,^)-box is ^‘^d.ii-h+c-h) _ 

• The index i measures the degeneracy of the defect of hyperbolicity. We have i = 0 


in the case of an initial ellipticity (Theorem 1.2), i = 1/2 in the case of a non semi-simple 
defect of hyperbolicity (Th eorem 1.3) and £ = 1 in the case of a semi-simple defect of 
hyperbolicity (Theorem 1.6). The instability is recorded in time 0(e| lne|)^/*'^'’"^) for initial 
frequencies 0(l/e). In particular, the higher the degree of degeneracy, the longer we need 
to wait in order to record the instability. 

• In the case i > 0, eigenvalues of the principal symbol are initially real (hyperbolicity). 
Instability occurs as (typically) a pair of eigenvalues branch out of the real axis at t = 0. 
The matrix Q should be understood as a change of basis, which includes a projection 
onto the space of bifurcating eigenvalues. The scalar fi corresponds to the real part of the 
bifurcating eigenvalues. Assumption 2.1 is formulated for the principal symbol evaluated 
along the bicharacteristics of /r. 

• In the non semi-simple case, the defect of hyperbolicity is typically not uniform in (x, ^). 
That is, if eigenvalues branch out of the real axis at initial time at the distinguished point 
(xq, ^o)) then the branching will typically occur for ulterior times t*(x, ^ > 0 for (x, ^) close 
to (xo,^o)- This is clearly seen in Lemma 5.1, under the assumptions of Theorem 1.3, and 


This is clearly seen in Lemma 5.1 
pictured on Figure 3. 

• The parameter q'*' corresponds to an upper rate of growth. In the elliptic case, is 
equal to the largest imaginary part in the initial spectrum, as seen in Section In the 
case of a smooth defect of hyperbolicity, q"'' = 9tA(0, x, ^), where A is a bifurcating 

eigenvalue, as seen in Section 


In the case i = 1/2, the block structure (2.7) derives from a reduction of the principal 


symbol to normal form; see Sections 5.1 and 5.2 in the proof of Theorem 1.3 


• In the case 1 = 1, the block structure (2.6) reveals a cancellation, seen on (6.4) in the 
proof of Theorem |1.6[ 


• The smoothly varying direction e (x) along which the lower bound (2.10) holds is not 
necessarily an eigenvector of A*; see the discussion in Section 2.3.3|and Lemma 5.10 


2.3.2. On verification of Assumption 2.1, We give in Theorems and |1.6| sufficient 

conditions, expressed in terms of the spectrum of A and the jet of the characteristic poly¬ 
nomial of A at t = 0 for Assumption |2.1| to hold. These sufficient conditions are satisfied 
in particular by Burgers, Van der Waals, and Klein-Gordon-Zakharov systems (Section]^. 
These conditions bear only on the coefficients of the system (the differential operator and 
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the source) and i?^(0), the initial datum of the reference solution. In particular, we may in 
practice verify these conditions without of course having any knowledge of for t > 0. 

2.3.3. On spectral conditions describing the transition from hyperbolicity to ellipticity. Con¬ 


ditions ( 1 . 8 ), ( 1 . 10 )-( 1 . 11 ) and ( 1 . 12 ) are all expressed in terms of the characteristic polyno¬ 


mial of A. Their generalizations in the form of conditions (2.10) and (2.11) are expressed in 


terms of the symbolic flow of A. Our point here is to explain why conditions bearing on the 
spectrum of A do not seem to be appropriate. The discussion below also highlights three 
difficulties in the analysis of the case £ = 1/2 : the lack of smoothness of the eigenvalues 
of the principal symbol, the lack of uniformity of the transition time (in the sense that the 
function t* does depend on (x,^)), and the lack of smoothness of the eigenvectors. 

A simple way to express the fact that an eigenvalue A of A branches out of the real axis 
at t = 0 is 

(2.13) A(0,x,^) G M, for all (x,.^) near (xo,^o)) with SAmA(0, xq, ■Co) > 0. 

But then by reality of the coefficients of A, eigenvalue branch out of M in pairs, so that 
( 0 ,xo,Co) is a branching point in the spectrum, and typically eigenvalues are not differen¬ 
tiable at a branching point, so that (2.13) is not nearly general enough. For instance, the 


eigenvalues of the principal symbol for the one-dimensional compressible Euler equations 

dtu + dxV = 0 , 
dtv + dxp{u) = 0 , 


(2.14) 


are 


A±(i,a;,0 = {u{t,x))Y^‘^. 

For a Van der Waals equation of state, for which there holds p' {u) < 0, for some tt G M, a 
transition from hyperbolicity to ellipticity occurs for data n( 0 , •) satisfying 

(2.15) p'(u(0,xo)) = 0, 9t(p'(u(0,xo))|^^g = -p"(M(0,Xo))9a:I^(0,Xo) < 0. 

The associated eigenvalues are 0(t^/^), in particular not time-differentiable at t = 0, so that 


condition (2.13) is not appropriate. 


A way around this difficulty is to consider the integral growth condition 


(2.16) 


A(r, X, i) dr = -f{t, x, 7(0, xq, ^0) > 0, 


for some i > 0 and some rate function 7 that is continuous in (t, x, ^) at ( 0 , xo,fo), and some 


local solution A of P = 0. Condition (2.16) may be verified with the Puiseux expansions 


of the eigenvalues at t = 0, such as in the Van der Waals example (for details in Puiseux 
expansions, see for instance chapter 2 in [ 8 ], or Proposition 4.2 in [25]). There are, however, 
at least two serious problems with (2.16). 


The first is that in (2.16), it is assumed that the loss of hyperbolicity occurs at t = 0 


over a whole neighborhood of (xo,^o)) which is typically not the case. Consider in this view 


the preparation condition (2.15) for the datum. From the second condition in (2.15), we 


find by application of the implicit function theorem that in the vicinity of (0, xq) the set 
{p\u) = 0} is the graph of a smooth map x —)• t*(x). The transition curve x —)> t*(x), defined 
locally in a neighborhood of xq, parameterizes the loss of hyperbolicity: for t < t*(x), there 
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holds p'{u) > 0 , implying QmX± = 0 , while for t > t^{x) there holds p'{u) < 0 , implying 
XsmX± 7 ^ 0. On the curve t = there holds p'{u) = 0, implying that the eigenvalues 

coalesce: A_ = A+. This means in particular that for (x,^) close to, and different from, 
we should not expect the loss of hyperbolicity to be instantaneous as in (2.16), but 
rather to happen at time and condition 2.16 should be replaced by 

(2.17) f ^mX{T,x,^)dT =-f{t,x,C){t-U{x,C)Y+^, 7 ( 0 , xq, Co) > 0, 

for some smooth time transition function t* > 0 , with t*(xo,Co) = 0 . 

The second issue with (2.16), still present in (2.17), is that while failure of time-differenta- 
bility of the eigenvalues is accounted for in ( |2.16 ), the associated lack of regularity of eigen¬ 
vectors is not. For instance, in the Van der Waals system (2.14), the eigenvectors of the 
principal symbol are e± := (1, ±(p'(u))^/^). In particular, under condition (2.15), the eigen¬ 
vectors e± are not time-differentiable at t = 0. It is then not clear how to convert conditions 
(2.16) and (2.17) into growth estimates for the corresponding system of partial differen¬ 
tial equations. Indeed, for instance in the simpler case of ordinary differential equations, 
spectral estimates such as (2.16) or (2.17) are typically converted into growth estimates for 
the solutions via projections onto spectral subspaces, an operation that requires smooth 
projections. 

We conclude this discussion by sketching a way around the issue of the lack of regularity 
of eigenvectors. Going back to the Van der Waals example, consider the ordinary differential 
equations 


J)®""’ 


= Id, 


parameterized by (x,^). Under condition (2.15), there holds 

p'{u{t, x)) = —a{x)t + 


for {t, x) close to (0, xq). Restricting for simplicity to the case p'{u{t, x)) = —t, we find that 
the entries (y, z) of a column of S satisfy the system of ordinary differential equations 

y' + iiz = 0 , z' - it^y = 0 , 
implying that y satisfies the Airy equation 


II j.c2 

y =t^ y, 


for which sharp lower and upper bounds are known. 

of the symbolic flow associated with the 


This motivates consideration, in Section 2.1 


principal symbol A. An important issue is then the conversion of growth conditions for the 
symbolic flow into estimates for the solutions to the system of partial differential equations. 
This is achieved via Theorem ID.31 


2.4. On the proof of Theorem |2.2| , We give here an informal account of key points in the 
proof of Theorem |2.2[ The proof is in three parts: (1) preparation steps which transform the 
equation into the prepared equation (3.36)-(3.37), (2) the use of a Duhamel representation 
formula, (3) lower and upper bounds. 
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(1) We introduce a spatial scale h and write perturbations equations about the reference 
solution (/>. We then block diagonalize the principal symbol (this is Q from the Assumption 
2.1), localize in space around the distinguished point xq, factor out the real part of the 
branching eigenvalues (this is symbol fj, from Assumption |2.1[ ) and change to a reference 
frame defined by the bicharacteristics of /i. Finally we operate a stiff localization in the 
elliptic domain D given by Assumption |2.1[ and rescale time. The key point in these 


preparation steps is to carefully account for the linear errors in the principal symbol, which 
take the form of commutators. The resulting principal symbol is a perturbation of the 


symbol A*(e,t,x,^) defined in Assumption 2.1 


( 2 ) Assumption | 2 . 1 | provides bounds for the flow of A*. As pointed out in Section 2.3 
these bounds bear on solutions to ordinary differential equations in finite dimension, in 
particular they are, at least theoretically, easier to verify than bounds bearing on spectra 
of differential operators. We use Theorem |D.3[ drawn from [2Tj and proved in Appendix [D| 
in order to convert these bounds into estimates for a solution to ( 1 . 1 ). 

Consider a pseudo-differential Cauchy problenij^ 


(2.18) 


dtu + op£(A.)ri = g, u{0) = uq £ L'^ 


where A is a symbol of order zero. Above, op£(A) denotes the e^-semi-classical quantization 
of symbol A, as defined in (3.1). Associated with the above Cauchy problem in infinite 


dimensions, consider the Cauchy problem is finite dimensions 


Theorem 


D.3 


5t5 + A5 = 0 , 5(r;r) = Id. 

asserts that if S{T]t) and its (x,^)-derivatives grow in time like exp( 7 A+^), 
with a rate 7 > 0 and a degeneracy index i > 0, then opg(5) furnishes a good approximation 
to a solution operator for 5^-|-op^(A), in time 0{\ lne|)^/(^+^^. That is, the solution of (2.18) 
is given by 


(2.19) 


u(t) ~ op£(5(OA))no/ oi>^{S{T;t))g{T)dT. 

Jo 


(3) The preparation steps (see (1), above) reduced our problem to a system of form (2.18). 
Via representation ( 2.19[ ), upper and lower bounds for the solution u to (2.18) are easily 
derived from the bounds of Assumption |2.1t and from postulated bound for the source g. In 
our proof, the source g comprises in particular nonlinear errors. Since we have no way of 
bounding solutions to 0 near (j) (the impossibility of controlling the growth of solutions 
with respect to the initial data being precisely what we endeavor to prove), we assume a 
priori bounds for the solution. The compared growth of opi,{S)uo and the Duhamel term 


from (2.19) eventually provide a contradiction. Note that the a priori bound (see (3.9) in 
Section 3.5) is particularly weak, since we allow for arbitrarily large losses of derivatives. 

We finally note that Garding’s inequality (see for instance Theorem 1.1.26 in m) asserts 
that nonnegativity of symbol A implies semi-positivity of operator op£(A). This is the 
classical tool for converting bounds for symbols into estimates for the associated equations. 


^^Notations and results pertaining to pseudo-differential calculus are recalled in Appendix B 
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It is shown in |24j how estimates derived from Garding’s inequality fail to be sharp in the 


non-self-adjoint case, as opposed to bounds based on Theorem D.3 

3. Proof of Theorem 12.21 


As discussed in Section 2.4, the proof decomposes into three parts: 


(1) Preparation steps which transform the original equation (1.1) into the prepared 
equation (3.36)-(3.37). This step covers Sections 3.1 to |3.10 


(2) The use of a Duhamel representation formula in which the solution to the prepared 
equation is expressed as the sum of the “free” solution (dehned as the action of an 
approximate solution operator on the datum) and a remainder - Sections 3.11 and 
[ 3 ^ 


(3) lower bounds for the free solution, and upper bounds for the remainder conclude the 
proof in a third step - Sections 3.13 to 3.15[ 

3.1. Initial perturbation. The goal is to prove ill-posedness, in the sense of Definition 
[IT} Parameters m,a,6,T are given, as in ( |1.6[ ), and we endeavor to disprove ( |1.7[ ). Define 

tpo{e,x) :=?fie(op^{Qe{0)~^)(^e^^''‘'^°^^''eeYx)y e > 0 , 

where {xo,fo) is the distinguished point in the cotangent space given in Assumption 2.1 
and 

• op£(-) denotes a pseudo-differential operator in e^-semi-classical quantization: 


(3.1) 


op^(a)u := (27r) 


-d 


X, e^^)D(^) df ,; 


Q£(0) = Q(0, xq + ^x, ^), with Q as in Assumption 


2.1 


• the vector e is as in Assumption 2.1 

• the spatial cut-off 0 G C, 
in .6(0,1/2). 


JOO 


has support included in B{0,S), and is such that 9 = 1 


Consider the following family of data, indexed by e > 0 : 


(3.2) 


u^(0, x) = 0(0, x) + 


X — Xq 

r-l — h 


where 0(0, x) is the datum for the background solution 0 (1.3), and K will be chosen large 
enough so that tt^(O) is a small perturbation of 0(0) in 77"* norm. 


Theorem 2.2 is a consequence of the following result for the family of initial-value prob¬ 
lems (l.l)-(3.2), indexed by e : 


Theorem 3.1. Given the parameters defined in (1.6), given a local solution 0 of (1.1) 
satisfying (1.3) with si large enough, under Assumption 2.1, if K is large enough: 

• either for all T and 6 with 0 < T < Tq, B{xo,5) C U there is no Cq > 0 such that 
for alio < e < Sq, the initial-value problem for (1.1) with the initial datum (3.2) has 
a solution in L°°([0, T], W^'°^{B{xo, 6)), 
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(3.3) 


or for some T and 6 with 0 < T < Tq, B{xo,6) C U, some Eq > 0, all 0 < e < eo; l^e 
initial-value problem for the system 0 with the initial datum ( |3.2[ ) has a solution 
in and there holds 

sup - - --- " --— = +CX) 

0<£<£o (0) ~ </’(0) ll//m([/) 

0<i<£'*T(£) 


where r(e) is defined in (2.3), so that, in particular, e"T(e) —)• 0 as e —>■ 0. 


Lemma 3.2. Theorem 3.1 implies Theorem 2.2 
Proof. There holds 


To 


X — Xq 
c-l—h 




X — Xq 


£ 


1-h 


where ip := opg((5(0, •,+ ')){&£), hence 

' X — Xq 


To 


-l-h 


—m-\-{l—h)d/2 


Hm(U) 


Let K > m — (1 — h)d/2. Then, 


a^(0, •) - (j){0, ■)\\h”^{u) 


^K—m-\-{l—h)d/2 


£—^0 


0 . 


Thus given a neighborhood U of i;/)(0) in and if e is small enough, then n^(0) lies in 

U. 


If for some £ small enough, the initial-value problem ( |l.l|)- (3.2 ) does not have a solution, 
then this means ill-posedness in the sense of Definition |1.1[ If there is a solution for any 
small e, then (3.3) disproves (1.7), since the sequence e^T(e) converges to 0, and again this 
means ill-posedness in the sense of Definition o □ 


3.2. The posited solution and its avatars. We assume that for some 0 < T < Tq, some 
0 < (5 with B{xo,5) G U, some Sq > 0, all 0 < e < Sq, the Cauchy problem for (1.1) with 
the initial datum (|3.2[)has a unique solution 


GL“([0,r],lTi’“(i?(xo,<5))). 


Our goal is then to prove (3.3). For future reference, we list here the successive avatars of 
the solution that we will use in this proof: 


u 

u* 

V 


perturbation 

spatial localization and projection 
convection 

stiff truncation and rescaling in time 


ii := {u‘^ — 4>){t,xo£^ ^x) 
:= op^{Qe){0u) 
u* := M*(0; t)v!’ 

V ■= {op^{x)u*){£^t), 


( ^ 

(M 

(3.27) 

(3.32) 
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3.3. Amplitude of the perturbation, limiting observation time and observation 
radius. The parameter K measures the size of the initial perturbation (3.2). We choose 
K to be large enough: 

(3.4) (2a — 1)K > 2am + (1 — a)(l — h)d, 

where m measures the loss of Sobolev regularity and a the loss of Holder continuity of the 
flow (as seen on target estimate (3.3)), and h = l/{l +£). 

The parameter T*, defined in (2.3), measures the final observation time in rescaled time 


frame. In the original time frame, the hnal observation time is (er*| lne|)^'^^^'*'^\ We choose 
T* to be large enough: 

(3.5) 7“(0,^o)7'* > ^, 

depending on K and the lower rate of growth 7 “ introduced in (2.9). 

The parameter 5 measures the radius of the observation ball H(0,(5) where the analysis 
takes place. (The radius is in the original spatial frame, and just 6 in the rescaled 

spatial frame associated with ft; see Sectionabove and ( |3.6| ).) If Theorem |3.1| holds for 
a given value of 5, then it holds for any smaller radius. In particular, we may assume that 
the given value of 5 is so small that the bounds of Assumption | 2 . 1 1 hold on |x| + |.^ — .^ol < 

In the final steps of our analysis (Sections 3.14| and |3.15[ ), we will further choose 5 to be 
small enough, depending on the growth functions 7 ^ introduced in Assumption |2.l| and T* 
(see condition (3.64) and the proof of Corollary 3.21). 

3.4. The perturbation equations. Our analysis is local in with 0 < t < e^T{e), 

\x — xo| < e^~^5 and — ^o| < where r(e) is defined in (2.3), and T* and 5 are defined 
in Section [TH 

The perturbation variable ii is defined in a rescaled spatial frame by 

(3.6) u{e,t^x) ■.= {u^ — with/i = 1/(1 + I"). 

The equation in ii is 

(3.7) dtii + £~^A{t, xo + £^~^x, £^dx)u + i?(e, t, x)u = F. 
where A is the 1-homogeneous principal symbol (|1.4|), 13 is order zero: 


B{£,t,x)u-.= 'Y^{duAj{t,Xo+£^ ^X,(j3e)u)dx^(t)e-duF{t,XQ+£^ '^X,4>e)u, 
j 

with notation (pe ■= Pit, xq + £^~^x). In d^ , the source F comprises nonlinear terms: 

(3.8) F = GQi£,t,x,u) ■ iu,u) + ^ Giji£,t,x,u) ■ {u,dxju), 

i<j<d 

where {u, v) —)• Go(e, t, x, ii) ■ (u, v) and {u, v) —>• G'ij(e, t, x, ii) ■ (u, v) are bilinear, defined as 

Goi£,t,x,u) := - [ (l-r)( V dlAjipe + ru)dxjpe - dlF{p^ + tu)) dr, 

•^0 i<j<d 

Gij{£,t,x,u) := — / duAjipi; + Tu) dr. 

Jo 














22 


NICOLAS LERNER, TOAN NGUYEN, AND BENJAMIN TEXIER 


We omitted above the arguments {t, xq + ^x) of d^Aj and d^F. In this proof, a pertur¬ 

bative analysis around (f> at {xo,Co), we will handle F as a small source, and B as a small 
perturbation of the principal symbol. 


3.5. A priori bound. The goal is to prove the instability estimate (3.3). We work by 
contradiction, as we assume that there exists C > 0 , such that for all t G [ 0 ,e^r(e)], there 
holds 


uniformly in (e,t), for 0 < f < e^T{e) = (eT*| lne|)^'^*'^''"^\ By choice of the initial datum 
( |3.2[ ), this implies (see the proof of Lemma |3.2[ ) 

(3.9) ||a'(() - 

with a possibly different constant C > 0, for all t < e^T{e). By definition of ii, this implies 

(3.10) \\u{t)\\whc^{B{o,S)) , fort<e^r(e), 

with notation 


(3.11) 


K' := a{K — m) — (1 — a)(l — h)d/2. 


By condition (3.4), there holds K' > K/2. 


3.6. Uniform remainders. The linear propagator in (3.7) will undergo many transfor¬ 
mations in this proof, through linear changes of variables corresponding to projections, 
localizations, conjugations, and so on. Every change of variable produces error terms. We 
will henceforth denote R^, for A: G Z, any bounded family Rfc(e,t) in in the sense that 


(3.12) 


sup ||Rfc(e,t)||fc,r < oo, 
0 <e<eo 
0<t<e^T{e) 


for r large enough, with notation || • H^,. for symbols introduced in (B.2) in Appendix]^ In 
the case /c = 0 , we say that a symbol belongs to Rq if either (3.12) holds or 


(3.13) 

By Proposition 


oo. 


B.l 


sup V sup |d“Rfc(e,t)|^i,jjd, < 
o<£<£o , 1 feM'* ^ 

o<t<£'^r(£) 

the corresponding operators op£(Rfc) are bounded 


L2 : 




L2 ^ 


\w\\e,s ■■= ||(1 + 


uniformly in 0 < e < Sq and 0 <t < e^T{e). We use above notation < introduced in (2.12). 
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3.7. Spatial localization and projection. The matrix-valued symbol Q{x, ^), introduced 
in Assumption |2.1[ is smooth, locally defined and invertible around (xq, ^o)- As explained in 
Appendix [C] we may extend smoothly Q into a globally defined symbol of order zero, which 
is globally invertible, with Q~^ G . We identify Q with its extension in the following, and 
let 

(3.14) u^{t,x) ■.= o\}^{Qe){0u), 


corresponding to a spatial localization followed by a micro-local change of basis. In (3.14), 
the function 9 = 9{x) is the spatial truncation introduced in Section 3.1, and we use notation 

(3.15) Qe{t,x,i) := Q(t, xq + (). 

Here opg(-) denotes a pseudo-differential operator in e^-semiclassical quantization, as in 
(3.1). Classical results on pseudo-differential calculus are gathered in Appendix]^ In par¬ 
ticular, opg(Q£) maps L?' to L^, uniformly in e, so that 

(3-16) ||u‘'||i2 < \\9u\\l^ < \\u\\(3(0,5))- 

We now deduce from the equation (3.7) in ii an equation in u'’, via the change of unknown 
(3.14). Here we note that the leading, first-order term in (3.7) is 

A{t, xo + e^~’"x, e’^dx) = op^{iAe), 


As := A(t,xo + £^ 


Thus the equation in is 


dtu^ +e ^op^{Qs)op^{iAs){9u) + op^{Qe){9Bu) - op^{{dtQ)e){9u) 
= oPsiQe)iOF) - ^ oi)^{Qs){Aj{4)s)udxj9). 


i<j<d 


At this point the goal is to express the terms in Oil above in the form of terms in u^, modulo 
small errors - that is, to approximately invert (3.14). This is done as follows. 

By composition of pseudo-differential operators with slow x-dependence (Proposition 


B.3), there holds 


Id = opJQ^ )oPe(Qe) + eopg(R- 


(3.17) 

where R_i is a uniform remainder in the sense of Section 3.6 With (3.17) we may thus 
express 9u in terms of u*' : 

(3.18) Oil = op^{Q~^)u^ -|- eop£(R_i)(0u). 

Using inductively (3.17), and composition of pseudo-differential operators (Proposition |B.2[ ), 
we obtain 

(3.19) Oil = op£(Ro)u*' -|- e”'op£(Ro)(0'u), 

for n as large as allowed by the regularity of (f>. 

By (3.18), the first-order term in the above equation in v!’ is 

oY>s{Qe)oY>e{'i^e){Ou) = op£(Q£)op£(fAe)op£((5“^)u‘' -b eop£(Q£)op£(iAe)op£(R_i)(6'M), 
implying, with Proposition B.3t 

OY)s{Qe)oVei}^e){Gu) = OY>s{iQe^eQe^)xt + £Op^(Ro){9u). 
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Besides, with ( |3.19[ ), we may write 

e^{oY>^{Qe){0Bu) - o^^{{dtQ)e){0u)) = 
where G Rq. From the above, the equation in li}’ appears as 

(3.20) dtv!’ + e~^OY)s{iQ£^£Qs^)'u^ + op£(R^)n*' = F*', 
where 

(3.21) B^ G Ro, := op,(Q,)(0F) - op,(Qe)(^,(<(>e)«5,,^0) + e”op,(Ro)(0n). 


3.8. Advected coordinates. Let M be the flow of opg(z;Ue) (or rather, as argued at the 
beginning of Section 3.7 of op^Ole), where /I is a globally defined symbol extending /j.) in 
the sense that 


dtM = op^{ifis)M, = Id, 

where the symbol jjL is introduced in Assumption |2.1[ and p,e is defined from p, by rescaling 


space as in (3.15). Let M* be the associated backwards flow, defined by 


drM* = —M*op(,(i^e), = Id. 


By hyperbolicity (reality and regularity of fi, and Proposition B.l), both M and M* map 
to L^, uniformly in e,t, for t < e^T{e). Egorov’s lemma (see for instance Theorem 4.7.8 
in m, or Theorem 8.1 in [21]) states that 

(3.22) M*M = Id+eop^(R_i), MM* = Id+eop^(R_i), 

where we recall that R_i is a generic notation for bounded symbols of order —1 (see Section 
3.6); in other words the equalities in (3.22) really mean that both M*M —Id and MM* —Id 
belongs to the class of operators of the form eop£(R_i). By Egorov’s lemma, given a G 5™, 
there also holds 

(3.23) M*op^{ae)M = op^{aei,) + eop£(Rm-i), 

where Og* denotes the symbol Og evaluated along the bicharacteristic flow, in the following 
sense: for any symbol q^we denote 

(3.24) K{t,x,^) 

where the (x*, ^*) is the bicharacteristic flow of /ig. What’s more, the remainder Rm-i above 
has expansion 

(3.25) op£(Rm_i) = op£(a*i) + eop,g(a*2) H-h e"'op,g(o*n) + e"+^op£(Rm-n-i), 

for n as large as allowed by the regularity of //, a and (j), where a*i G 5"™“* has support 


included in the support of a*. Identities (3.23) and (3.25) also hold if M and M* are inter¬ 


changed, with backwards bicharacteristics replacing forward bicharacteristics. By (3.22), 
there holds 

M*op,(Ro) = M*op,(Ro)(MM* + eop,(R„i)), 


^'^Except for symbol A, as seen on definition on T* in (2.4l. 
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implying, by ( 3.23[ ): 

M*op£(Ro) = op£(Ro*)M* + e(M*opj(R_i) + op£(R_i)M*), 
and reasoning inductively we arrive at 

(3.26) M*op£(Ro) = opg(Ro)M* + e''op£(Ro), 

where n is as large as allowed by the regularity of (f). The advected variable is defined as 

(3.27) M* := M*(0;^)u^ 

There holds 

dtu^ = M*{dt - op^{ii^e))u\ 

and dfU^ is given by equation (3.20) in u^. Using ( |3.22 ), we find 

M*Op^{iQeiAe - fJ,e)Q7^)v!' = Op^{{iQe{As - fJ,e)Q~^)^)u^ 

+ eop£(Ro)u* - £M*op^{QsA^Q~^)op^{'R^i)u\ 


in 


with notation (3.24). Thus, with (3.26), 

M*op^{iQeiAe - = Op^{{iQe{Ae - 

+ eop£(Ro)u* + e”opj(Ro)u^ 

Besides, in view of (3.26), the order-zero term in ( |3.20 ) contributes to the equation i 
Ui, the terms 

M*(0; t)op^{B^)v!’ = op£(R*)m* e”op£(Ro)u^ R* G Rq. 

The equation in u* thus appears as 

(3.28) dtUi, + £~^op^{{iQe{Ae - ^e)(57^)*)u* -h op£(R*)m* = F*, 
where 

(3.29) R* G Ro, R* := + e’^op,(Ro)u^ 

The symbol {Qe{A^ — /re)Qj^)* is symbol Qe{As — fJ,e)Qs^ evaluated along the bicharacter¬ 
istics of /ig, as defined in ( 3.24[ ). 

3.9. Frequency, space, and time truncation functions. We introduce frequency cut¬ 
offs XoiXo,Xo^ spatial cut-offs 9\^,9o,9o, and temporal cut-offs 'ijjQ,'tpo, 'ipQ. All are smooth 
and take values in [0,1]. Given two cut-offs and 1 /^ 2 ) 

(3.30) V’l ^ ^2 means (1 — V’ 2 )iAi = 0. 

Equivalently, -< ^1^2 when ■02 = 1 on the support of i/'i. 

The supports of the frequency cut-offs Xo:Xo,Xo all assumed to be included in the 
ball {|^| < 5}. All three are identically equal to 1 on a neighborhood of .^ = 0. There holds 
Xo^Xo^ Xo- 

The supports of the spatial cut-offs 9q,9o, 9q are all assumed to be included in {|x| < h}. 
All three are identically equal to 1 on a neighborhood of x = 0. There holds 0 q -< 9o -< 9o. 


We assume in addition 9q < 9, where 9 is the spatial cut-off of Section 3.7 
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Figure 6. The truncation function i/^o- 


The temporal cut-offs i/lQji/jo) '*/’o are nondecreasing, supported in {t > —(5}, and identi¬ 
cally equal to one in a neighborhood of {t > 0}. In particular, i/jq = 1 on {t > —(5/3}. There 
holds i/o ^ ^0 ^ V’o- The truncation i/o is pictured on Figure 6. 

Associated with these cut-offs, define 

(3.31) := xo Oo{x)tpo{t - f*(e,x,0), 

and dehne similarly x^ and x ia terms of XQ,0g,V’o and respectively. In (3.31), 

the transition function f* is defined in (2.1), which we reproduce here: 

t*(e,x,0 = with 6»* > 0, 6i*(0,^o) = 0, Va;,^6»*(0, ^o) = 0. 


Recall that the value of C is fixed in Assumption |2.1[ depending on £ : there holds C = 0 if 
.£ = 0 or = 1, and ( = 1/3 li £ = 1/2. 


Lemma 3.3. The support of x is a neighborhood of the elliptic domain V (2.2), and there 
holds 

\dxd^x\ < 

where C is introduced in Section 2.1.1\ such that, in particular, f < h. The same holds of 
course for the other truncations x~^d x, which satisfy X^ ^ X ^X- 


Proof. For {t, x, £,) to belong to the support of y, there needs to hold simultaneously |C~?o| < 
|x| < 5, and t*(e, x,^) — 5 < t. This defines a neighborhood of V (precisely, of the 
projection of V onto the (t,x,^) domain), as defined in (2.2). 

We may now assume t* to be not identically zero, otherwise x is not stiff. Then (see 
Section 2.1.1) £ = 1/2, h = 2/3, f = 1/3. By the Faa di Bruno formula 


di 


- t*(e,x. 




(k) 


01 -h/3fc =0 


(t-t*) n 


-hfi0: 


9^0*), 


where are positive constants. We note that 5|*0*(e,x,^) = 0{e‘i) if |ci:| = 1 and f 

belongs to the support of xo; by assumption on 0*, while 9/^0* = 0(1) if |q:| > 2. 

Consider the case of a decomposition of /3 in a sum of ftj of length one. Based on the 
above formula and the bound on the corresponding bound is 
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If /3 is decomposed into /3i + /32 + • • • + /3fc, with |/3i| = 2 and |/3j| = 1 for j > 2, then 
k = IP] — 1. The corresponding bound is < ^-l/^K^-C) as soon as C < h/2, 

which holds true. It is now easy to verify that the decomposition of |/3| into sums of multi¬ 
indices of length one corresponds to the worst possible loss in powers of e. 

We turn to x-derivatives of — t*). By assumption, 0* is a function of 
and e^~^{x,^). Thus x-derivatives bring in either powers of since h = 2/3, 

or powers of £-^+i-^+C = since = h/2 = 1/3. 

Thus ~ ^*)| ^ if = 1/2. Considering finally the full truncation 

function y, we observe that the term in xo contributes the exact same loss per ^-derivative, 
and the spatial truncation Oq contributes no loss. □ 

Corollary 3.4. The operator op£(x) maps to uniformly in e, and so do 

op£(x^) and op£(x). 


Proof. Since x is compactly supported in x, we may use pointwise bounds for 5"x a-iid 
bound (B.5) in Proposition B.l The result then follows from Lemma 3.3 □ 


3.10. Localization in the elliptic zone and rescaling in time. We define 
(3.32) V := op^{x{t)){u^{eH )), 

meaning that we first rescale time in n* and then apply op£(x) evaluated at t, where x is 
defined just below (3.31). We now derive an equation in v, based on equation (3.28) in ?x*. 

Consider first the leading, first-order term in (3.28). When evaluated at e^t, its symbol 
is precisely (2.4) the rescaled and advected symbol for which Assumption 2.1 holds: 

H'iQsiAs - iae)Q7^)*) L x, f). 

Thus 


OPe (x) (opg {{{iQe (Ae - ^ )*) j [eH) = op^ (x)op£ (iA*) (u* {eH) ). 

Similarly, denoting B := Bi,{E^t), where is the order-zero correction to the leading 
symbol which appears in equation ( |3.28D , there holds 

op£(x)(op£(B*)M*)(e^t) = opg(x)op^(B)(u*(e^t)). 

We now introduce a commutator: 

-h EB){Ui,{E^t)) = op^{iA^ -)- eB)v + r(rt*(e'"t)), 

where 


(3.33) f := [op^{x),op^{iA^ + eB)]. 

By definition of v, 

dtv = op^{dtx){u^{EH)) + E^op^{x){{dtu^){EH)). 

Together with equation ( |3.28| ) and the above, this implies 

dtv + E^~^op^{iA^ + eB)v = op^{dtx){Ui,{E^t)) - £'*“^f(u,*(e^t)) -h £^op^(x)F*. 

We will handle the right-hand side as a remainder. The following Lemma shows that we 
may introduce a truncation function in the above principal symbol. 
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Lemma 3.5. For any bounded family P{s,t) G S^, there holds, for x defined in (3.31) and 
V defined in (3.32).• 

(3.34) 0Ve{P)v = 0Ve{xP)v + [op^ (P), opg(x)]i^ + e” (op^ (Ro)ix*) (e^t), 

where Rq is a uniform remainder in the sense of Section \i.Q\ 

Proof. By definition of x and x, and Proposition B.2[ there holds 

OPe(x) = opfixx) = OPe(x)oPe(x) +e”'^OPe(P„/(X,x)), 

and the remainder satisfies 

l|op£(Pn'(x,x))llL2^||.||,,_„, < \\d'^'x\\o,C{d)\\d^'x\\o,C{d)- 


We use here norms || • \\m,r for pseudo-differential symbols of order m, as defined in (B.2). 
By Lemma |3.3[ there holds 


Xllo,C(d) 


< p-{P+C(d))C 


and 


9” xllo,c(d) < 


Since Q < h, there holds e'^'^ in'+ 2 C{d)){h Q ^ j£ jg chosen large enough. 

Thus 

(3.35) op£(x) = op£(x)opg(x) + e''op^(Ro), 

where Rq is bounded for t < T{e). This implies 

o^fiP)v = o^fiP)o^fix)u^{£H) = opg(P)opJx)^' -t-£''opjRo)(u*(e'"t)). 

Now 

opfiP)opfix) = opfixP) + [op£(P),oP£(x)], 


and (3.34) is proved, with a symbol Rq which is a uniform remainder in the sense of Section 
3.6, meaning that we rescale in time the remainder which appears in (3.35). □ 


Applying Lemma 3.5 to P = iA* -|- eB, we derive the final form of the equation satisfied 
by V : 


(3.36) 


dtv + e^ ^opfixiiA^ + £B))v = £^g, 


where the source term g is defined in terms of the remainder F^, from (|3.28|)-(|3.29|): 

(3.37) 

g = e“i(ru - f (u*(e'*t))) -h £-’^opfidtx){u*{^'"t)) + oPe(x)(P*(e^i)) + e”(op£(Ro)«*)(e''t), 


with T is defined just like T (3.33), but with x in place of x- 

The derivation of (3.36)-(3.37) ends the first step of the proof of Theorem 2.2 Our goal 
is now to show a growth in time for the solution to (3.36) over the interval [0,T(e)], where 

T(e) = (T*| 11161)^^^^"*"^^ In this view, we will first derive an integral representation formula 
for V. 
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3.11. An integral representation formula. At this point we use the theory developped 
in Appendix [D) Theorem D.3 gives the integral representation formula for the solution v to 
( |3.36[ ) issued from u(0) : 

(3.38) u = opg(S(0;t)))u(0) f op£(i:(r; t))(Id+eop£(Ro))(^(r) + eRou(0))dr, 


where Rq are uniform remainders, as defined in Section [3. 6 t and the approximate solution 
operator op£(S(s;t)) is defined by 

(3.39) Yj 

where qo is large enougtj^ In (3.39), the leading term 5o is defined for 0 < r < t < T{e) 
and all (x, ^) G x by 

(3.40) 9fSo + + eR)(t,x,05'o = 0, S'o(r;T) = Id, 

and the 5^, for (7 > 1 , are correctors, defined inductively as the solutions of 

(3.41) 

dtSq + e^~^x{iA^ + £B){t,x,()Sq + e^~^ ^ (xiiA^ + eB)%^Sq2, S’g(r;r) = 0. 

91+92 =9 
0<gi 


with notation crijlnO '2 ■= (—i)”(n!) ^ X]|a|=n From (3.40) and (3.41) we deduce 

the representation, for q > 1, 

(3.42) Sq{T;t) =f S’o(r';t) ^ (x(iA* + eR))(r', x, 0tlgi-S'q2(r; r') dr'. 

91+92=9 

0<9i 


In order to be able to exploit representation (3.38), we need to check that Assumption D.l 


under which Theorem D.3| holds, is satisfied. This is the object of the forthcoming Section. 

3.12. Bound on the solution operator. Recall that S, the symbolic flow of is 

defined in (2.8), and is assumed to satisfy the bounds of Assumption |2.1[ The upper bound 
(2.11) in Assumption 2.1 is assumed to hold for S in domain D defined in 


:= {(r;t,x,0, < t < t < T{e), |x| < d, |^-^o|<d£^}, 

where the transition time is defined in (2.1) and the final observation time T(e) is defined 
in (Q. 

The goal in this Section is to show that the symbolic flow Sq, which is defined as the 
solution of (3.40), and the correctors Sg defined as the solutions to (3.41|), satisfy the bounds 


of Assumption [D . 1 This will allow to use the representation Theorem |D. 3 [ and will justify 
representation |3.38 This will also give a bound for the norm of the approximate solution 
operator op£(S) defined in (3.39). 

We are looking for bounds for the correctors Sq, and their derivatives. Consider repre¬ 
sentation (3.42). Disregarding (x, ^)-derivatives, we see on (3.42) that Sg appears as a time 


11 


Depending on the final observation time T* (2.31, and the growth function 7 ; see Appendix [d| and 


in particular the proof of Lemma 


D.2 
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integral of a product 5o(x(*^* + £B))Sq^, with q 2 < q. We may recursively use (3.42) at 

rank (72; and by induction Sq{T‘,t) appears as a time-integral of form 

(3.43) 

/ So{Ti;t){x{iA-, + eB)){n)SoiT2-,Ti){x{iA^ + eB)){T2)--- 

J T<Tl<"-<Tq<t 

■ ■ ■ So{Tn-, rn-l)(x(i^* eB)){Tn)So{T-, Tn)dTl . . . dTn, 

in which there are q occurences of iA^, + sB and O' -|- 1 occurences of Sq. Note again that in 

(3.43) we overlook (x, ^)-derivatives. From there, it appears that we need to 

• derive bounds on Sq and its (x, .^)-derivatives; these will be deduced from the bounds 
of S postulated in Assumption |2.H 

• derive bounds for products of 5o with -)- eB); here the block structure as¬ 

sumption 2.6-(2.7) from Assumption 2.1 will come in. 

3.12.1. Product bounds for S. In a first step, we prove bounds for products of symbols 
involving the symbolic flow S (2.8) and the rescaled and advected principal symbol A* 
(2.4). We denote for a,l3 £ and 0 < r < t : 

SaA'^'A) := e^~^S{T-,t)d^d^A^{T), 

and products, for {aiAi) £ and 0 < r < < t„_i < ti < t : 

(3.44) Sji(t, Ti, ... ,Tn]t) := (ti; t)S 02 , 1^2 ('^2! Ul) • • • iTn] Jn— i)5’(t; Tn). 

Lemma 3.6. Under Assumption \2.1\ there holds 

1 


|S„(r,ri, ... ,Tn\t,x,i)\ 


rC 


1 


e +(r;t,x,0, 


for all n > 1, all a*, jdi £ all 0 < t < Tn < ■ ■ ■ < ti < t with (r; t, x, ^) G V, uniformly 
in £. By < we mean here entry-wise inequalities modulo constants, for each block o/S„, as 
described below (2.12). 


Proof. If all blocks of A* satisfy (2.6), then ^ = 0, and the stated bound simply follows 
from the multiplicative nature of the growth function, namely identity 

(3.45) e^+(ri;t)e.^+(r;ri) =e.y+{T;t). 


Suppose then that a block A*j satisfies (2.7), and consider the associated component Sj of 
the symbolic flow (2.8). Bound (2.11) from Assumption 2.1 and the cancellation observed 
in product 


(3.46) 


1 e 

,1-h 


h-l 


1 e 

,1-h ^ 


h-l 


= 2 


1 £ 
A-h ^ 


h-l 


imply, omitting the index j, 

(3.47) |5(ri;t)9“5^A*(ri)5(r;Ti)| < 


£ 


1 £‘ 
1—h 


h-l 




The result follows from (3.45)-(3.46)-(3.47) by a straightforward induction. 


□ 
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3.12.2. Product bounds for S. 


Next we show that spatial derivatives of S, and products 

We denote for a, /3, 


involving d^S, satisfy the upper bound (2.11) from Assumption 2.1 
G 0 < T < t : 


Sa,i3{T-,t) :=e^ ^d^S{T]t)d^df A^{t), 
and products, for a, G (3i G and 0 < r < < r„_i < ri < t : 

Sn(T, Ti, ... ,Tn'it) := {t\ ; t)Sa2,fii2 (”^2) Tl) • • • Tn—l)d^ S{T',Tn). 

Lemma 3.7. Under Assumption \2.1l there holds the bound in domain T) : 

1 


|5“5| + |S, 


< 

r\j 


gC 1 1-7-’ 

for each block of the block-diagonal matrices S and S^. The precise meaning of< is described 
below (2.12). 

Above, and often below, the time and space-frequency arguments are omitted. In partic¬ 
ular, the “interior” temporal arguments of S^, namely r„,... ,ri, are omitted. It is implicit 
that the r* are constrained only by r < < Tn-i < ■ ■ ■ < ti < t, and that the indices ai,j3i 

and n > 1 are arbitrary. 

Proof. We first prove by induction on |a| that d^S enjoys the representation 
(3.48) <9“5 = ^ f Sn, 

n<|Q| 

where S„ is defined in ( 3.44[ ). By (3.48), we mean precisely 
d^S{T-,t)= a/ 

-1 ^ ^1 I r 


(3.49) 


Sn(T, , . . . , Tfi, t) dri . . . dTn, 


l<n<|a| 


' T<Tl<-"<T„<t 


with constants Cn independent of (e,t, x,^). In the following, whenever products, such as 
S„, are integrated in time, the integration variables are the “interior” variables, as described 
just above this proof. In order to prove (3.48) for |q;| = 1, we apply 9“ to the equation (2.8) 
in S : 

dtdfS + ie^-^A^dfS = -ie^-^{dfA^)S. 

This implies the representation 

df.S{T\t) =—ie’^~^ f 5(r';t)9"A*(r')S'(r;r')(ir', |a!| = 1, 


which takes the form (3.48). For greater values of |q;|, we have similarly 

d-S = -ie^-^ Y. f Sd^^A^d^^S, 

OLi-\-a.2=cx. ^ 

0<|qi| 


and the induction step is straightforward. From (3.48) the bound on dfS follows by the 
bound on Sr, in Lemma 


3.6 


and the multiplicative nature of Time integrals only 
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contribute powers of |lne|, which are invisible in < estimates. Finally, from (3.48) we 
deduce (using the same notational convention as in (3.48)) 

(3.50) S„ = > ' I Sr 


s„= 

n'<|ai|H-|-|a„|+n+|a| ' 


and Lemma |3.6| applies again. 


□ 


3.12.3. Bounds for the symbolic flow Sq. We bound here Sq, the solution to (3.40). While 
S is the flow of the symbol Sq is the flow of + sB), where x is a stiff 

truncation and i? is a bounded symbol of order zero. We prove here that the upper bound 
(2.11) in Assumption |2 .1 1 is stable under the perturbations induced by x and B, in the sense 
that So and its spatial derivatives satisfy the same upper bound as S. 

In a first step, we consider the solution to 

(3.51) dtSy^ + e^~^xA*S^ = 0, S^{T;T) = ld. 

Associated with S^, define products involving 9“5^ just like S„ was defined as a 
product involving d°S, but with x^* in place of A*, explicitly: 


S^^a, 0 {r-,t) := e t)a“2a[(xA*(r)), a = (ai,a 2 )GN ,/3 G N , 

and products, for cij G fli G and 0 < r < < r„_i < ri < t : 

Sx,n('E, Ti, . . . , Tn'i t) := (ui) 0'S'x,a2,/92 Ui) • • • Sy^^cii„,l3n (Ai! Tn-l)dx Ai)- 

Corollary 3.8. The solution to ( |3.51| ) enjoys the bounds 

1 


(3.52) \d^S^{T;t,x,^)\ 
and 

(3.53) <£-1^1^ 
for any 0 < t < t < T{s), any (x,^). 


r-C 

1 




e +(r;f,x,0. 


Just like the bounds of Lemma 3.7, the bounds of Lemma 3.8 are understood entry-wise, 
for each block of the block-diagonal matrices and 

While the bound of Assumption |2.1| was stated over domain P, the bound of Corollary 


3.9 holds for any {x,ff), for 0 < r < t < T{e), with T{e) defined in (2.3). This comes from 


the truncation x in (3.51). 
Proof. There are hve cases: 


If T < f < f* — (5, then x = 0 on [r,f] (see the definition of x in (3.31)), implying 


S^ = ld. 


• If r < f* — J < t, then by property of the flow and the previous case 

S’x(r;f) = S^{T-,t^)S^{U-,t) = 5^(4; f), 
and we are reduced to the case t^, < t < t < T{s). 
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If < T < i < r(e), then x = 1 on [r; t], and = S hy uniqueness. 

If t* — 5 < r < t : comparing the equation in S (2.8) with the equation in 


(3.51), we find the representation 
S’x(r;t) = 5(r;t) 

and, applying 5“ to both sides above, 
(3.54) 


- x)(r')A*(r')5'x(r;r')dT', 


dxSxir; t) = d^S{T-, t)-e^ ^ ^ t)d^^ ((1 - x)(r')A(r'))9""5x(^; ^0 dr', 

where the sum is over all ai + a 2 + <^3 = ol. Since we are only interested in upper bounds. 


we omitted multinomial constants > 0 in (3.54). We factor out the expected growth by 
letting := e“+i 

solves, omitting the summation sign, 

rt 


(5v, S'^ := e.iS'. Then, by property (|3.45|) of the growth function, 


(3.55) 5“5^(r;t) = a“5 (r;t)-e' 


A-l 


ar5'(r';t)9r((l-x)(r')A(r'))5“3 5^(r;T')dr'. 


Now given a matrix M = 


mil rni2 
m 2 i m 22 


we let M := 


mil e^mi2 
e“‘^m2i 17122 


. There holds 


identity, analogous to (3.46): 

(3.56) Ml M 2 = M 1 M 2 . 

In particular, if M^ and M o are bounded in £, then M 1 M 2 is bounded in e. Thus from 
(3.55)-(3.56) we deduce 

The key is now that e^“^cl“((l — x)4*)9 is uniformly bounded in e. Indeed, by Lemma 

S spatial derivatives of the truncation x are uniformly bounded. By the block conditions 
-(2.7) and definition of just above (|3.56), the matrix £^~^A^ is uniformly bounded 
in e. Besides, 5 is uniformly bounded, by Lemma 


3.7 


Thus we obtain the bound 


|a“4(r;t)|<C(l + I |5“35^(r;r')|dr'), 

with 1 0 : 3 1 < I a I, for some C > 0 , implying by Gronwall and a straightforward induction 
|5“5^(r; t)\ < U - 5 <t <t <U, 

which is good enough since t — t < 5. Going back to we find bound (3.52). 

• The same arguments apply in the remaining case t* — (5<r<t*<f. 

At this point (3.52) is proved and we turn to (3.53). First consider products involving no 
spatial derivatives of Sy^, which we denote for consistency with notation (3.44). Here 
we note that the proof of Lemma 3.6 uses only the upper bound ( 2 . 11 ) for S, via cancellation 


(3.46). We may thus repeat the proof of Lemma 3.6 and derive a bound for Sy^n- The only 


difference is that, while A* is uniformly bounded in e, the truncation function x is stiff in 
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as seen on (3.31) and definition of in (2.1), and reflected on Lemma (3.3). This gives 
bound (3.53) with in place of 

Finally, spatial derivatives are well behaved, in the sense that they are bounded without 
loss in e. Thus from the bound in we derive bound (3.53) exactly as in the proof of 
Lemma 3.7, via a representation formula identical to (3.50), with and in place of 
Sn and Sn. □ 


Corollary 3.9. The solution Sq to (3.40) enjoys the hound 

|5"5'o(r;t,x,0l < e"‘^e.^+(r; t, x, 0, 

for any 0 < t < t < T{s), any 


Proof. From equations (3.51) in Sy and (3.40) in 5o, we deduce the representation 


(3.57) 


9"5o(r;t) =5“5x(r;t) + e^ / d^"S^{T'-,t)d^^{xB{T'))d^^So{T-,T') dr', 


with an implicit summation in 01 + 02+03 = O; and implicit multinomial constants 
Cat > 0. We factor out the expected growth be fore applying Gronwall’s Lemma, as we let 
5 q := e“+S'o and := From Corollary 3.8 we know that < £~^. Besides, 

from Lemma 3.3 and the fact that B G Rq, the symbol d°^{xB) is uniformly bounded. 
Thus from (3.57), we deduce 

\dxSo{T-,t)\<£-‘^ + £^-‘^j |5“3S’^(r;r')|dr', I 03 I < |o|. 

We may now conclude by Gronwall’s Lemma and a straightforward induction, since r(e) 
grows at most logarithmically and h — (^ > 0 . □ 

3.12.4. Prod uct bounds for Sq. The next step is to prove product bounds for O^Sq, as in 
Lemma 3.7 In this view, we let for a = (ai, 02 ) G /3 G N'^, 0 < r < t : 

and products, for a* G /3j G N'’* and 0 < r < < r„_i < ti < t : 

(3.58) So,n('r, TI, . . . ,Tn]t) := S’OjOi,/?! (ui; t) 80 ^ 02 ,f }2 Ul) • • • So^an,l3n {'^n'l Xn—l)d^ Soir^Tn). 

Corollary 3.10. There holds the bound 


(3.59) 


< p-C(l+l/3|), 


for any 0 < r < < Tn-i < • < Ti < t < r(e), any (x,^), any at, fit, any n > 1, with 

l/5| = EilAI- 


we 


Proof. First observe that, in the case n = 1, if we were to directly use Corollary 3.9 
would find the upper bound 

£^-Vr^o(r';t)<9r A(r')9“^5o(r;T')| < 

which is not good enough since h — 1 — (j <0ifC = l — Hence a need for more than the 


bound on 9"5o from Corollary 3.9 
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Second, note that the loss in (3.59) is coming from ^-derivatives applied to X) exactly as 


in the proof of Corollary 3.8 


^x,k + 


cB 


|a| > 0, 


We are going to use the representations 

(3.60) d^So = 

for any m G N, with notation 


and 


Sf,o,fc := E^\d^^S^d'i^{xB)) • • • 






In (3.60) we use compact notation as in (3.48). In particular, time arguments are implicit 

and form “chains” in the sense that in products {S^dx^xB){S^dx^xB), the time arguments 
of the first term are (T;ri) and those of the second term are (Ti;r 2 ). The integrals in 

(3.60) are time integrals bearing on the “interior” variables of and S^; see the explicit 
definition of S in (3.44) and see how compact notation in (3.48) is expanded into explicit 
notation in (3.49). 


By Assumption 2.1 and Corollary 3.9 there holds 


(3.61) 




< .Kh-O.-Cr 


Representation (3.60) is proved by using recursively (3.57) in itself n — 1 times. 
From (3.60), we find So,i = e^~^d^^Sod^^d^{xAi,)d°'^So to be a sum of 9 terms: 
The term is bounded with Corollary 

For terms II = 


3.8 


{xA-^)d^^S^, we use Corollary 


3.8 


for the rightmost product, 
of form e^~^dx'°^^ Sy-d^^d^ {xA^,)d^^ Sy., and then bound separately the remaining k products 
of form £^dx*Sdx' (xB)- This gives |II| < and we use ( < h. 

• Term III = may not be handled by the same argument as 

II, since the rightmost product here has form Sod'^'^d^{xA^)d^^Sy_. Here we use 

(3.61) with k = m large, the bound \d'^'^d^{xAic)\ < , and the bound of Corollary 


3.8 


for There occurs a loss of but this is compensated by an appropriate choice of 

m. Precisely, we find |III| < e™-(^“f)“2C-|/3|Ce^+ and then choose m large enough, depending 
on /i, so that m{h — C) ~ C ^ 0- 

• Term IV = £^~^d^^S^d^^d^{xA^,)S^i. is symmetric to II. Here we isolate the leftmost 
product (xA*)(9“^ 5^. 

for 


In term V = e 


— ^S^j^d^^d^{xA-i,)S^j^,, we use the cancellation of Corollary 


3.f 


the term in the middle Sy^d^^d^{xA^)dx^ Sy^. The remaining terms contain k + 
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k' occurences of S', and there is a prefactor sK^+k') _ Thus we obtain the bound |V| < 
^(fc+fc')(/x-C)e-C(i+l/9|)e^+. 

• The remaining terms all involve at least m factors, hence, by (3.61), have an 
prefactor. We handle these terms just like III above. 

From the above, we conclude that there holds 

|So,i| < 

where |/3| is the number of ^-derivatives in So,i. 


The general case n > 2 is handled in exactly the same way. Via representations (3.60), 
products involving O^Sq, as the statement of the Corollary, are expanded into sums of 
products involving 5“ S^, and remainders which involve products with a large number of 
terms, and So terms. These remainders are handled as term III above, using h < Q, 
hence as small as needed for m large. We are then left with products involving only 

X-^-ki xB and spatial derivatives of S*^. For these, we use Corollary 3.7 as we did above in 
the treatment of terms I and II (using the case n > 2 in Corollary 3.8, while we used only 
n = 1 in the above treatment of terms I and II). □ 

The final bound in these preparation steps involves products of Sq with as 

already met in (3.43). We let for a = (ai, 02 ) G /3 G 0 < t <t : 

:= -he5)(r)), 

and products, for ctj G /3j G N'^ and 0 < r < < r„_i < ti <t : 

(3.62) S^„(r,ri,... ,r„;t) := t)5o,„2,/32(r2; n) • • • T„_i)a"5o(r; r^). 

Corollary 3.11. There holds the bound 

loR I < .-C(l+l/3|)p ^ 

for any 0 < r < < Tn-i < • < ti < t < T{s), any (x,^), any ai,f3i, any n > 1, with 

l/5| = EJAI- 


Proof. Developing, we find to be a product of terms of form So,n' (3.58) with terms of 
form E^dfd^xB. By Corollary 3.10, there holds |So,n| ^ where |/3| is the total 

number of ^ derivatives that appear in So,n- Besides, \e^dfd^{xB)\ < , The result 

then follows by C < h. □ 

3.12.5. Bounds on the correctors Sq. We are now ready to prove bounds on spatial deriva¬ 
tives of the correctors Sq introduced in (3.41): 

Corollary 3.12. There holds the bounds for 0 < t < t < T{e), for any {x,f) : 

, for0 <q<qo. 


Proof. For (7 = 0, Corollary |3.9| gives the desired bound. For q = 1, the corrector Si admits 
representation (3.42), so that d^Si appears as a sum of terms 

rt 


^h-1 


d^^Soir'-, (x(zA + eB))(r')a“3+-5o(r; r') dr', 
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where |a| = 1, and ai + 02 + as = a. Corollary 3.10 applies and gives the desired bound. 
Consider now representation 3.42 for Sg, any g > 2. In this representation, d°Sg appears as 
a sum of terms 


(3.63) 


-h-l 




where \k\ + q' = q, |fi;| >0, and ai + a 2 + as = a. We may recursively use representation 
(3.42) in (3.63), so that d^Sg appears as a time integral of terms (3.62), with exactly 
q derivatives bearing on the ^ variables. The result then follows from Corollary 3.11[ □ 

3.12.6. Bounds on the approximate solution operator op£(S). We arrive at a bound for the 
action of the approximate solution operator opg(S) defined in (3.39): 

Corollary 3.13. There holds the bound, for 0 < r < t < T(e) : 


|oPe(S(r; t))w\\L2^Bixo,S)) ~ ^ ^ 


\w 


where := max 7 (x,£). 

|x|<<5 

l«-Col<<5 

Proof. Let 01 be a spatial cut-off that is identically equal to one on a neighborhood of 
B{xo,d). Then 

\\oPsi^iTlt))M\L2(B(xo,5)) - 

Now by Proposition B.l[ 

l|oPe(0iS(T;t))u;||i2(Rd) < ^ supJ5"(0iS(t; t, •, 0) |ii(Rd) 

|a|<d+l 


By Corollary 3.10, there holds for all 0 < g < qo, for |^ — ^o| < d : 

e‘'^|0i5“Sq(r;t,-,OlLi(Rd) < sup e +{T]t,x,(,). 

\x\+\^-^o\<s 

The result then follows from (^ < h and the pointwise bound 

e.^+(r;t, x,^) < exp( 7 +(t^+^ - ), for £ > 0. 


□ 


3.12.7. Lower bound for Sq. We verify that the lower bound (2.10) in Assumption 2.1 
stable by perturbation: 


IS 


Lemma 3.14. For 6 and e small enough, the flow Sq satisfies lower bound (2.10) from 
Assumption 2.1, that is 

(0; T(e), X, ^ 0 ) < | Sq{0;T{£), x, ^o)e{x) \ , 
for |x| < 6 and e(x) as in Assumption\2.1\ 
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Proof. First observe that t*(e, 0, ^o) = 0, so that the equation in S and coincide over the 
time interval [5,T(e)], at (x,^) = (0,^o)- A simple perturbation argument, similar to the 
arguments developed in detail above, then imply that 5^(0; T(e), 0, ^o) satisfies the lower 
bound (2.10). Next we use representation (3.57) from the proof of Corollary |3.9[ with a = 0, 
which implies 

|*S'oe| > \S^e \ t)||.B(r')||5'o(r;r')| dr'. 


As argued above, we may use lower bound (2.10) for besides. Corollaries 3.8 and 3.9 
provide upper bounds for Sq and S^. These yield the lower bound 

|S’o(0;r(e),x,^o)e| > e"^e^-(0;r(e),x,^o) - e^"‘’e^+(0; r(e), x, ^o)- 
Now at (r;t, x,^) = (0; r(e), x, ^o)) there holds for d and e small enough 

- niin 7"(-,^o)) 

\ ur <(5 


e.v+e^ f < exp 


■'7-^7- 


The constant 


do := max 7 +(-, ^o) - min 7 (', 6 ) 
l?-«o|<<5 

is small for small d, by continuity of 7 ^ and the fact that 7 ^( 0 , ^ 0 ) = 7 ~ ( 0 , fo)- In particular, 
we may choose d, depending on T* introduced in (2.3) and satisfying condition (3.5), so that 

(3.64) h — ( — T*do > 0. 

Thus for |x| < d and e small enough, there holds 

1 -e''-^+e-i(0;r(e),x,eo) > 1 - > 1/2, 

and the result follows. □ 


3.12.8. Bounds on d^d^Sq. We finally give bounds on (x, ^)-derivatives of So and correctors 
Sq. This ends the verification of Assumption D.l 


Lemma 3.15. There holds bounds, for 0 < t < t < T{s), for any (x,^) : 


(3.65) 




I < p-C(i+l/3|+g), 


Proof. From representation 

d^So{T-,t) =^ / 5o(T';t)9^^(x(iA*+ eS))(r')9^^5o(T;r')dr' 

hi+h 2 =h 

l/3i|>0 


and Corollary 3.11 we find by induction bound (3.65) for q = 0. Besides, applying 
the equation in Sq, we find the representation 


d^d^Sq = -e 


h-1 


So ixiiA, + eB))d^^dl^Sq + ((x(fA* + ) ), 


with an implicit summation over 01 + 02 = o, /3i + /32 = /?, (?i + (72 = O' with |oi| > 0 , 
\fii\ > 0, (71 > 0. We use recursively the above representation, and find dfd^Sq to be a sum 
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of terms with a total number of ^ derivatives equal to |/3| + q. It now suffices to apply 
Corollary |3.11[ □ 


3.12.9. Conclusion. Before moving on to the third and last part of the proof of Theorem 

verify Assumption 


3.15 


D.l 


Thus at this 


2.2, we recapitulate our arguments so far. 

The bound and Lemma 

point the integral representation (3.38) is justified, via Theorem D.3 We reproduce here 
equation (3.38): 

V = op£(S(0;t)))u(0) [ op^{T.{T-,t)){ld+£op^(Ro)){g{T) + eop^(Ro)v{0))dT. 

Jo 

This ends the second part of the proof. Next we will prove a lower bound for v, in norm 
L'^{B{0,S)). In this view, we will bound from above the time-integrated term in (3.38), and 
bound from below the “free” solution op£(S(0; f))u(0). The bound from above for the time- 
integrated term rests on Corollary 3.13 bounding the action of op£(S), and on a bound for 
the source g, which is the object of the forthcoming Section 3.13[ The bound from below 


for the free solution is a consequence of the above Lemma 3.14 


3.13. Bound on the source term. Our goal in this Section is to give an upper bound 
for the source term g from equation (3.36). The source g is defined in (3.37), which we 
reproduce here: 

g = e-i(ri; - f -h £-'^op^{dtx){u^{£^t)) + op£(x)(T;(e^t)) -b e”(opg(Ro)«*)(e''t). 

We let 

90 ■■= op£(x)(^*(e''i)) + e’"(op£(Ro)ii*)(e^i), 

and 

g^ := e-^Tu - f -b £-'^op^{dtx)M£H)), 

so that g = go + Ox- We first consider g^. 

Lemma 3.16. There holds the bound 


\\oPeix'')9o{t)\\L^ < ||ft(e''t)|lwi.oo(B(0,5)) +e''||ft(e^t)||L2(B(0,5))> 

uniformly int G [0,r(e)], for n as large as allowed by the regularity of cf, where the possibly 
stiff truncation x^ is defined just below (3.31), in particular such that x^ A X) with notation 
introduced in (3.30). 

Proof. There are three types of terms in go : (a) a commutator issued from the spatial cut¬ 
off 6 , (b) nonlinear terms, (c) remainders of size 0(e"'). The term T* is defined in (3.29), 
in terms of defined in (3.21), itself defined in terms of T (3.8). We reproduce here their 
definitions: 


F = 


Go{£,t,x,u) ■ {it, it) -b 




F^ = op^{Qe){9F) - ^ op^{Qe){Aj{(j)s)udr,.e) -b e”Ro(6'tt), 
F* = M*F^ + e"op,(Ro)n^ 
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(a) Commutator term: this is the term in dx^O in F''’. 

This term is rendered small by the left action of op£(x^)- Indeed, commutators that arise 
from a localization step depend on derivatives p', where p is the localization symbol. If we 
further localize with p*', such that is identically equal to 1 on the support of p, then from 
(I — p)p^ = 0 we deduce p^p' = 0, and the associated commutator is arbitrarily small. This 
is made precise below. 

Consider first the term 


OY>Ax)^*OVe{Qe){^Me)udxF). 

We start by approximately commuting op£(x*') and M*, by use of (3.23) (with M in place 
of M* and conversely): in view of Lemma 3.3, there holds 

(3.66) op£(x^)M* = M*opj(x(*)) + e^"^opg(Ro), 

where x\i,) denotes evaluation of along the backward characteristics of ps- Note that we 
use here definition (3.13) for Ro. In particular, the characteristics depend on time through 
e^t. Expanding in powers of e^t, for t < T{e), we find terms that are all supported in 

the support of x^ up to a remainder 0(e”), with n as large as allowed by the regularity of 
cj). In conjunction with Proposition |B.2t this implies 

M*op^{x\i,))oPeiQe) = op^iOoKo) + e"op£(Ro), 

where 9o is defined in Section 3.10, in particular such that (1 — 9)6q = 0. Now 

op^{9(iIlo){Aj{(l)s)udxj9) = op^{9odxj9Ilo){Aj{4>s)u) + [op£(6»oRo), dxj9]{Aj{(l)^)u), 

and the first term above is identically zero, since 9Qdxj9 = 0. For the second term, we use 
Proposition B.2: up to 0(e”), the operator involves product of 9q and its derivatives with 
derivatives of dx 9. These products are identically zero. Thus 


\\M*Op^{x\^^)op^{Qe){Aj{cl)e)u)\\L2 




Now Egorov’s lemma ( |3.23[ ), as used in (3.66), may be expanded up to arbitrary order, as 
in (3.25). The supports of the symbols that appear in this expansion share the property 


that we used for x\^y We obtain 


\\0Peix'’)M*0p^{Q^){Aj{(pe)udx^9)\\L2 < e"||lj||L2(B(o,5))- 
(b) The nonlinear terms are local in ii and ’VxU, so that 

9{x)Go{e, t, X, ii) ■ {ii, ii) = 9{x)Go{e, t, x, 9'^u) ■ {9'^u, 9'^u) 


where 9^9^ (see notation (3.30)), with supp0* C A similar identity holds for Gi. 

Thus 

\\9Go{u) ■ {u,u)\\l^ < \\9Gq{9^u) ■ {9h,9h)\\L2 

<G{\9^u\L^)\9^u\L^\\9h\\L^. 

Since 0** has support in there holds ||0^n||i2 < ||?j||L2(5(o,(5))) same in L°°. By 


the a priori bound (3.10), there holds in particular 1111^00(^(0^5)) < 1, since K' > 0. Thus 
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A similar argument for Gi yields an upper bound that involves |^t|^yl,oo(^(o,(5))• We conclude 
that 

hVe{X^)^*OVe{Qe){0F)\\L2 < ||h||^i.„o(B(o,5))- 

(c) The remainders of form e”opj(Ro)(0'it) and e"'op£(Ro)u^ in g con tribu te e”||0t6||2,2 
to the estimate for opg(x^)( 7 , by definition of uniform remainders (Section |3.6[ ), and bound 
(3.16) in v}’. Same for the remainder in u*, since ||u*|| 2,2 < ||u^|| 2 , 2 , by property of M*. We 

□ 


use in addition Corollary |3.4| for all these remainders. 

Corollary 3.17. There holds, for any P G S^, for the source term go defined just above 
the statement of Lemma 3.16|, the bound 


.2K' 


(3.67) l|op£(x‘')op^(P)go||L2 < \\P\\o,C{K'), 

uniformly int G [0,T(e)], for some C{K') > 0, and where || • ||o,r is the norm in defined 
in (B.2). The constant K' is defined in ( 3.11[ ) in terms of K,a,m and d. 

Proof. By Lemma |3.16| and Proposition |B.1[ there holds 


|opg(P)op^(x^ 


Il2 < 


\P 


lo,C(d)(ll«llL2(S(0,5)) + £”11^11^2(5(0,5))), 


and with (3.10) we obtain the upper bound (3.67) by taking n = K '. There remains to handle 

modulo terms which are 


B.2 


the commutator [op£(x^), opg(P)]. Here we use Proposition 
0(e"’), the symbol of the commutator is a sum of terms of the form P^ := ~ 

d^X^d^P)- Sine x^ P X (using notation introduced in (3.30)), there holds Pa = xPa- Thus, 
by Proposition B.2 again 

(3.68) op^{Pa)g = op^{Pa)op^{x)g + e^^^°‘^+^'>op^{Qa)g, 


where the leading terms in Qa have the same form as Pa- For the first term in (3.68) above, 
we use Corollary 3.16, as we may since X < X- We use also the fact that ||oP(,(Pq,)|| 52 _^ 2.2 is 
bounded uniformly in e, by Lemma 3.3 and bound (B.5) in Proposition B.l[ For the leading 
terms in Qa, we use inductively (3.68), and arrive at (3.67). The Taylor expansion in the 
composition of operators needs to be carried out up to order 0{K'), hence a dependence in 
||P||o,c(A'') iu the upper bound. □ 

Lemma 3.18. There holds the bound, for any P G S^, for g^ defined just above Lemma 

[3T61- 

I|oP£(x^)oP£(P)5xIIl 2 < e^'+”||P||o,n, 
uniformly in t G [0,r(e)], for n as large as allowed by regularity of cj. 


Proof. First consider the case P = Id. S ince derivatives of fjo vanish identically on the 
support of i/’q, there holds, by Proposition 


B.2 


\0Pe{x'')0Pe{dtX)\\L^^L^ ^ x‘'llo,C(d) l|5x9fX||o,C(d)- 


By Lemma 3.3 this implies 

l|op,(x')op,(aix)||L2^L2 < e-^e-(-+2C(d))c^ 


and the above is arbitrarily small if n is large enough, since (j < h. The same argument 
holds for the other two terms in g.^. In the general case of P G 5*^, we may reason as in the 
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proof of Corollary 3.17 that is by first spelling out the composition up to a 


large order, then using the above. Remainders are small by condition Q < h. 
The above argument prove the bound 

(3-69) l|op£(i/^^)op^(P)5^||i,2 < e'^{\\u^{eH)\\L2 + \\v\\l2 ). 


We finally use bound (B.4) in Proposition B.l in order to control ||u||^ 2 (]gd), and obtain, via 
Lemma |3.3| 

lkllL2(]Rd) 

This loss is irrevelant since in (3.69) the integer n can be chosen to be very large. We 
conclude by bound ||u* IIl 2 < ||ij^||l 2 and bound ( |3.16[ ) for □ 

3.14. Lower bound for the free part of the solution. First we describe the time 
transition function and the truncation ipQ for frequencies close to ■ 

Lemma 3.19. There holds the bound, for |x| < <5 and |^| < 

0 < U{e, X, ^0 + < Cd, 

for some C > 0 independent of S. In particular, for such (x, ^), there holds 
'0o( - U{e,x,f,o + e’^C)) = 1, V'o(r(e) - U{£,x,Co + e'^0) = 
and same for 'ipo and V’o- 

Proof. We may assume i = 1/2 here. There holds 

t*(e, X, ^0 + 

= Co + 

Jo 

and, by assumption on 0* (2.1), for |x| < <5 and |C| < there holds 

\£-^e^{£^-^x,Co)\ < ^ ^^ 2 ^ 

and 

|d^6»*(e^“^x,Co + e^uOl < C'(e^“'*|x| +e^|C|) < Ce‘^ 6 , 

for some C > 0 which does not depend on 6 . This proves the bound on t*. In particular, for 
such (x,C), there holds — t * > —6/9, and T(e) — t * > —6/9, implying the result for i/q, by 
definition of i/q (see for instance Figure 6). □ 

Based on the above result for t* and ipo, we identify the leading term in the datum for v 
defined in (3.32): 

Corollary 3.20. There holds 

(3.70) ||u(0) - E^P^<°/^eo{x)e{x) ||^, < 

where 9q is the spatial cut-off introduced in Section 


3.9 
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3.2 


Proof. At this point the reader may hnd useful to jump back to Section 
Id, the datum for v is 

^;(0) = e'^op£(x(0))opg(Q£(0))f6i3f?e fopg(Q£(0)"^) 


Since = 


We may comm ute o Pr(Qs) and 9, since this produces an error that is in thanks 

to Proposition B.2 Then in the datum we handle separately the oscillations in and 


the oscillations in —fQ/e^. There holds 

e^op,(g,(O))op,(Q,(O)-i)(e'(->«°/^'*0e) = + 0(e^+^), 

where O(-) denotes a control in L^. Thus the oscillation in f,o/£^ contributes to u(0) the 


term 


£^op,(y(O))(e-«°/^'^02g-') = e^e-«°/^^p,(x(0,x,eo + •))(^'e) 

modulo terms that are 0{£^^^). By dehnition of x ia Section 


3.9 


In the above we may use Lemma 


x(0, X, ^0 + = Xo(e^ ‘^O^o{x)i’o{-Ui£, x, 

since Xo(^^~‘’0 is non-zero only if |^| < 6 e^~^. Thus 


3.19 


x(0,x,^o+ = Xo(e'* = 0o(a;)^0(0- 

We may expand the above up to an arbitrary power of e^“^. The remainder is a stiff 
symbol in f , hen ce we will lose in evaluating its operator norm (in accordance with 

Proposition B.l), but such a loss is irrevelant if the prefactor is large enough. Also, 

the Taylor expansion brings out a large prefactor, implying that the norm of the 
remainder depends on the high Sobolev norm ||6*oe 11^?^, but Oq G (7“ and e is assumed to 
be smooth, hence this derivative loss is irrevelant as well. We thus obtain 

By choice of 9o, there holds 9o -< 9. Thus 9o9‘^ = 9o, and we obtained the leading term 
£^9o{x)e (x) as claimed in (3.70). 

It remains to prove that the oscillation in has a small contribution to the datum. 

The leading term in the datum associated with this oscillation has the form 

+ ■)Pe)0, 9 := W, 

where G S^, uniformly in e. The key is then that op£(x(— .^o + which is smooth 

and supported in P(0, <5) (because 9o is supported in B{0,6)) is pointwise bounded by 
l'^MLi(|c|>c/P)’ ^ ~ ^l^ol ~ e^(5/2 > 0. This norm is arbitrarily small, since T9 

belongs to the Schwartz class. Spatial derivatives are handled in the same way. □ 


Corollary 3.21. For V defined in ( |3.32[ ), there holds for e and 6 small enough the lower 
bound 


||op£(x^(r(e)))op^(Po(0;r(e)))u(0)||^2(^(Q_5)) > ‘^exp(7 r(e)^+^), 

where 7 “ := min| 3,|<5 7 “(x, ^ 0 ) o.nd y*' is introduced in Section 


3.9 
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Proof. According to Corollary 3.20[ the datum is 

u( 0 , x) = e^v{0, x) + u( 0 , x) = e^^'^°^^'^9o{x)e{x). 

There holds 

(3.71) op,(5o(0;r(e)))i}(0)(x) = r(e), x, ^ 0)^0 e» + eVo, 


with 


Fo ;= ^ op, 

|a|=l 


{d^So) ( 0 ;r(e),x,^o+ 'r(-))d'r (5“(6»oe)) (x). 


We now apply op,(x^(T'(e))) to (3.71). The leading term is 

e“-^o/^'*op,(x^T(e), X, ^0 + •))(So(0; T(e),eo)0oe). 

By definition of Section 


3.10 


x‘'(T(e),x,Co + = Xo(e^ ‘^CWoix)f’o{T{£) - U{e,x,Co + s'^0)- 

We use Lemma |3.19t as we did in the proof of Corollary |3.20[ there holds 

x‘'(T(£), X, ^0 + = 6»^(x) ^1 + ^ ^ dr 

The leading term now appears as 

Uo = e“'«°/^''5o(O;r(e),x,eo)0^(x)e(x), 

since 9q -< 9q. And since 0g = 1 in a neighborhood of 0, we may use Lemma 
states that for s and 6 small enough and |x| < 5, there holds 


3.14 


which 


£ ^e, 


^-(0;r(e),x,^o) < I 5’o(0;r(e),x,^o)e(x)|. 


Consider the lower growth function e..^-, as defined in (2.9), at (x,^) = (x,^o)- It involves 
t*(e,x,^o)- By Lemma 3.19, 0 < t*(e,x,^o) < Cd. Thus there holds the lower bound 


r-C 


exp 


(7 T{£)^+^)<£ ^e.^-( 0 ;r(e),x,^o), 


uniformly in |x| < 6 , with 7 “ defined in the statement of this Corollary. Hence a lower 
bound for the L‘^{B{0,6)) norm of the leading term of the free solution as desired. 

There remains to bound from above the terms we overlooked so far. The first invo lves 
the r emainder in the datum, which is 0 {£^~^^~^) in norm. With Corollaries |3.4| and 
3 .I 3 I the action of op,(5o) on this remainder is controlled by £^~^£^~^ exp ( 7 '’'T(e)^"’'^). 


The other terms are e^+^op,(x*^(T(e)))Vb and ‘’Wo, with notation 

Wo := ^ op, (9|'xo(e"‘’'^(-)) (5“(5*0(0; T(e), •,^o)0oe)) (x). 


|a|=l 


By Corollary 3.4 Lemma |3.15| and Proposition |B.1[ there holds 


.K+h 


|oPe(x‘'(P(e)))Vb||L 2 (Rd) < ^ exp ( 7 +T(e)^+^). 
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For Wo, we observe that op£(5|*Xo) ^ Fourier multiplier, bounded in —)• norm by 

the maximum of its symbol, equal to a constant independent of e. We may then bound 9“S'o 
with Corollary 3.9 Thus e'^+^“^Wo is controlled in just like e^+^opg(x^(r(e)))Vb- 


Summing up, we obtained a lower bound of form 

exp (7-r(e)^+^) (l - exp ((7+ - 7-)r(e)^+^)). 

By definition of T(e) in ( |2.3[ ), 

exp ((7+ - 7-)r(£)i+^) = e/^-C-(T+-7-)T*_ 

Given T*, since h — C > 0, we may choose <5 small enough so that the difference 'y~^ — 7“ 
is so small that h — — 7“)r* is strictly positive. The result follows is e is small 

enough. □ 


3.15. Endgame. We apply op^{x^{T{e))) to the l eft of the representation formula (3.38) 
for V at t = T{e), with defined in Section 3.10, and prove that the contribution of the 


initial datum dominates the time-integrated Duhamel term. This eventually provides a 
contradiction to the assumed a priori bound (3.9), and concludes the proof. 

Based on (3.38), we find 

\\op,ix\T{e)))viTiemL2^Bm) > I - (H + HI)- 

The leading term is 

I = oPe(x‘’(T(e)))op^(5o(0;T(e)))u(0). 

This term is bounded from below in <5)) norm by Corollary 

I|I||l2(s(o, 5)) ^ e-^“^exp(7-r(e)^+i). 

The error term in the contribution of the datum is 


3.21 


We control II by Corollary 


II := 5'“’op,(x‘(T(E)))op,(S,(0;r(£)))l<0). 

(action of op£(x*') in T^), Lemma 


3.4 


3.15 


their derivatives) and bound (B.5) from Proposition |B.l This gives 

|II||l2(s( 0,5)) < e-^+'^'^exp (7+r(e)^+^). 


(bounds for Sg and 


The Duhamel term is 
III = e 


oPeix\T{e)))op^{j:{T-T{£))){ld+eop^{'Ro)){g + eop^{'Ro)v{0))dt'. 


10 


We bound III with Corollary |3.17 and Lemma 3.18[ in which we choose n = K' : 

I|III||l 2 (b( 0,5)) < exp (7+r(e)^+^) (e^^' -b . 


Above, we used the bounds of Corollary 3.12|in o rder to control the norm ||L!||o^(g(j^ /p w hich 
appears in the upper bound of Corollary 3.17 Since 2K' > K hy condition (3.4), we 


obtained 


|op,(x^)u(T(e))||i.(B(o,5)) > e^-^exp exp {j+T{eY+^). 
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We may now conclude as in the proof of Corollary 3.21 Precisely, rewriting the lower bound 
as 

l|oPe(x‘')^^(r(e)) 11 ^ 2 ( 5 ( 0 , 5 )) > e^-‘^ exp ( 7 “r(e)^+^) (l - , 

and choosing 6 small enough so that h — ( 7 "*“ — 'y~) T^ > 0 (recall that 7 ^ is the local 
maximum of the rate function 7 from Assumption | 2 .l| and 7 “ is the local minimum of the 
lower rate function 7”), we hnd for e small enough the lower bound 

\\oVeix^)v{Tie))\\L2(^B{0,5)) > C\ 

for some C > 0 independent of e, where | lne|* is some powers of | lne|. By choice of T* in 


(3.5), there holds K — — ^ T* < 0 if 5 is small enough. Hence a lower bound which blows 


up as e —J 0, contradicting the a priori bound (3.10). 


4. Proof of Theorem 11.21 initial ellipticity 


We are going to verify that, under the assumption of initial ellipticity (1.8), Assumption 


2.1 holds with parameters 


i = 0, h = 1, C = 0 ) and = 0 , = 0 . 


Then Theorem 1.2 is a consequence of Theorem 2.2 


4.1. Block decomposition. Let Ao,Ai,...,Ap be the spectrum of A at (0,xo,^o)- The 
ellipticity condition (1.8) states that one of the Xj is non real. By reality of A, complex 
eigenvalues come in conjugate pairs. In particular, at least one of Xj has strictly positive 
imaginary part. We may assume Q'mAo > 0 and 


(4.1) $5771 An > max SjmA,-. 

i<i<p 

Let m be the algebraic multiplicity of Aq in the spectrum, and Eo{t,x,^) the generalized 
eigenspace associated with the family Ao,i(t, x, ^),..., Ao,m(L of (possibly non distinct) 
eigenvalues of A which coalesce at (0,xo,'Co) with value Aq, that is Ao,j(0, xq, Co) = ^o- Let 
i?i(t, x,C) be the direct sum of the other generalized eigenspaces. Let Qq be the projector 
onto Eq and parallel to Ei. The Ao,j may not be smooth, but are continuous (see for instance 
Proposition 1.1 in [25]), and Qo is smooth (see for instance | 8 |, or Proposition 2.1 in |25j). 
and determines a smooth change of basis Q{t,x,^) in which A is block diagonal: 

(4.2) QAQ-‘=:('"f J . 

We focus on the block A(o) associated with eigenvalues Xqj. The symbolic flow accordingly 
splits into S'(o), where S'(o) solves 


Ot^{0) + ^H(O) [et, Xo + X, 4)5(0) = U, 5(0) (r; r) = id . 

By (4.1) and a repeat of the arguments below, the component 5(i) of the symbolic flow is 
seen to grow not as fast as 5(o) near ( 0 ,Co)- 
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4.2. Reduction to upper triangular form at the distinguished point. Let now P be 

a constant (independent of t, x,^) change of basis to upper triangular form of A(o)( 0 ) 
and Qfj^ be the diagonal matrix 

= diag(l,..., 

The parameter /r will be chosen small enough below. There holds 

QnP^{o){^^xo,Co)P~^Q'il^ = AoId+/iJ, 

where J is upper triangular, bounded in /r, with zeros on the diagonal. By a Taylor expan¬ 
sion of A(^Q-^{et,xo + x,^) in {t,x,^), we observe that there holds 


i^(o)(et,xo + x,C) = P ^Q^^{iXold+inJ + B{e,t,x,i))PQ^, 
where the Taylor remainder B has form 

B = £tBi{e, t, X, 0 + (x, ^ - ^o) • B 2 {£, t, x, 0, 
with Bi and B 2 bounded, uniformly in e, in domain 

(4.3) |x| < (5, t<r*|lne|. 

Let 

(4.4) S{T;t) := exp (j,{t - T)\o^PQf,S(^o){T-,t). 

Then S solves 


(4.5) dtS + (i^iJ + B)S = 0, S{t]t) = PQ^. 

4.3. Bounds for the symbolic flow. Consider the Hermitian matrix 

3f?e {ifiJ + B):=^ {{i^J + B) + {ifiJ + B)*). 


Its eigenvalues A are semisimple, and vanish at {fa,t,x,^) = (0,0,0, .^ 0)5 hence satisfy, for 
{t,x,^) in domain ( |4.3[ ), the bound (see for instance Corollary 3.4 in [25]): 

|A| < co(^ + er*| lne| |x| + - Col) =: 7 o(e, A^,LC), 


for some cq > 0 independent of e, fi, t, x, C, for {t, x, C) hr (4.3) and (e, small enough. Thus 


(4.6) 


—70 Id < iRe {ifiJ + B) < 70 Id . 


Let 


S± := exp ( ± (t - T) 7 o(e,/r,t,x,C))S. 


From (4.5) we deduce, for any fixed vector e G 


-'Si.e{dtS±e,S±e)crn + + B) ± Jo) S±e, S±e) = 0. 


By (4.6), this implies, for e unitary, 

(4.7) < |5(r;t)e| < 

Back to 5 ( 0 ) (4-4), we now have 

|5(0)(r;t,x,C)| < 
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We now choose // = | lne| and let 


7 + ;= 3=771 Ao + Co(|x| + 1C - CoD- 
We obtained, for s small enough, 

|5(o)(r;^)|<|lnere(*-")^^ 

corresponding to the upper bound (|2.11 ). Finally, since |(PQ^)“^e| > ci/i^“™'|e| for some 
Cl > 0 and all e, we deduce from (4.7), the lower bound 

\S{o){T-,t)e\> , 7 " = SmAo - co(|x| + 1C - Col), 

corresponding to the lower bound (2.10). 

5. Proof of Theorem II.3t non semi-simple defect of hyperbolicity 


It suffices to verify that, under the assumptions of Theorem 1.3, Assumption 2.1 holds, 
with parameters 

7 = 1/2, /i = 2/3, C = l/3. 

Then Theorem 1.3 appears as a consequence of Theorem |2.2[ We may assume initial 
hyperbolicity (1.9), since otherwise Theorem 1.2 applies. 

5.1. The branching eigenvalues. Let ojq := (xq, Ao,Co)- By assumption, there holds 

(5.1) P(0,a;o) = 0, dxP{0,u}o) = 0, dtPiO,uJo) / 0, 5|P(wo) / 0. 

By the second and fourth conditions in ( |5.1[ ) and the implicit function theorem, there 
exists a smooth function //*, with /7*(0, xo,Co) = Aq, such that d\P = 0 is equivalent to 
A = x,C), for (t,a:,C) close to (0,xo,Co)- 

By the first three conditions in (5.1) and the implicit function theorem, there is a smooth 
r*, with r*(a:o,Co) = 0, such that P(lt*) = 0 is locally equivalent to t = r*(x,C)- 

We now use the above implicitly defined functions /i* and r* to describe the spectrum of 
A near (t, x, C, A) = (0, wq) . 


Lemma 5.1. In a neighborhood of (0, wq), there holds P = 0 if and only if 

(5.2) (A - A7(x,C))^ = -{t - r*(x,C))e(Lx,C, A), 

where := ;^*(r*(x, C), C) c smooth and satisfies e(0, xq, Co,-^o) > 0- 


Proof. Given t close to 0 and (x, C, A) close to wq, we expand the characteristic polynomial: 

P{t, X, C, A) = P(r*(x, C), x, C, 1 i*(t*(x, C)), X, 0) 

(5.3) +{t- r*(x,C))ei(t,x,C) + (A - n{x,f)fe 2 it,x,f,,X) 

= {t- T*(x,C))ei(t,x,C) + (A - fi{x,f)fe 2 it,x,f,,X), 

since P(r*, •,/7*(r*)) = 0, with 

ei(t,x,C) := / (<97P)((1-T)r*(x,C)+ rt,x,/7(x,C))dT, 

Jo 

e2(t,x,C, A) := [ (1-r)(5fP)(t,x,C, (1-r)/7(x,C)+ rA)fir. 

Jo 
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We let e := 6162 . Then e(0, xo,^o, Aq) > 0, as a consequence of the definition of ei and 62 


and condition (1.11). The result follows from (5.3) 


□ 


Equation (5.2) describes a pair of eigenvalues branching att = r*(x, out of the real axis, 


with imaginary parts growing like {t — The time curve t = r*(x,^) is the boundary 

between the hyperbolic region t < in which the eigenvalues are real, and the elliptic 

region t > r*(x,^) in which the eigenvalues are not real and where we expect to record 
an exponential growth for the symbolic flow. In the introduction, Figure 3 pictures the 
hyperbolic and elliptic zones in the {t,x,^) domain near (0,a:o,^o)- 

We note that under the assumptions of Theorem |1.3t the above defined time transition 
function r* satisfies 

(5.4) T* > 0, T(((xo,^o) = 0, r"(xo,Co)>0, 

where T((xo,Co) is the differential at (xo,Co) and r*(xo,^o) the Hessian. Indeed, if the first 
condition in (5.4) were violated, then hyperbolicity would not hold at t = 0, contradicting 


(1.9). Thus 0 is a global mininum for r*, and (5.4) ensues. 


Remark 5.2. Note that under assumption d\PdtP < 0, eigenvalues stay real for small 


t > 0, by Lemma 5.1 


5.2. Change of basis. By Lemma 5.1, the eigenvalues A±(t,x,^) of A{t,x,f,) — /i(x,^)Id 
satisfy 

(5.5) A±(t,x,0^ = -{t - n{x,^))e{t,x,C,p{x,C) + X±{t,x,^)), 

where e is smooth in all its arguments, and e{0,xo,Co, hi^o^^o)) > 0, with //(xo,^o) = Aq. 
From (5.5) and continuity of e, we deduce the fact that A_ and A-|- are purely imaginary, 
hence A+ + A_ =0, since the matrix A — //Id has real coefficients. From ( |5.5[ ), we also 
deduce the bound, for some C > 0, locally around (0,xo,^o)) 

|A±| < C\t-U\^/^, 

which we may plug back in ( |5.5| ) and deduce, by regularity of e, 

(5.6) A± = ±i{t - T*)^/^e(t, X, //(x, 0)^^^ + 0{t - r*). 

We now reduce A to canonical form: 

Lemma 5.3. There exists a smooth change of basis Q such that locally around (0,xo,^o), 


Q{t,x,^){A{t,x,f) - //(x,0Id)Q(t,x,0 ^ = 


A 


( 0 ) 


0 


0 A 


( 1 ) 


where Hfoi = ( /_, x ^ \'^/2 > o,nd A(i\ G 2 )x{N 2) smooth. In 

\ -{t-T^)eo + o{t-T^yp oy’ 

the bottom left entry o/H/g), the funetion r* is the time transition function introduced just 


above Lemma 5.1, and we denote eo(t, x,^) = e{t,x,f,,p{t,x,f,)), with e as in Lemma 5.1 


The block associated with Ag is size two, since the multiplicity of Ag is equal to two (|1.10|). 

Bo 


0 


0 A 


( 1 ) 


Proof. We may smoothly block diagonalize A — p, for instance as described in Section 4.1 

he 

Thus for some smooth Q, there holds Q{A — p)Q~^ = 


where Rg is the 2x2 
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matrix Bn = 


d 


ail ai2 


, with smooth entries aij. The spectrum of Bo is {A_, A+}, where 

021 022 J 

X± satisfy ( |5.6[ ). Since, as noted above, there holds A_ = — A+, the trace of Bq is zero, that 
is 022 = —On- Besides, there holds 021012 / 0 at (0, xo,Co)- Indeed, if 021012(0 , xq, Co) = 0, 
the spectrum of Bq would be smooth in time, contradicting (5.2). Without loss of generality, 

\ 


we assume 021(0, xq, Co) 7^ 0. Then Qq = 


0 O21 

1 “O21 Oil 


is a smooth change of basis such 


that QqBqQq ^ = 


0 1 
A 0 


The bottom left entry of Bq is equal to the product A_A+ of 


its eigenvalues. By (5.6), we find that A_A+ = —{t — T*)eo + 0{t — r*), and the result holds 

Qo 


with Q = 


0 

0 Id{^Ar-2 


Q. 


□ 


5.3. The symbolic flow. Our goal is to prove the bounds of Assumption 2.1 for solution 
S to the ordinary differential equation 

dtS{T] t) + re“^/^A*(t)5(r; t) = 0, 5(r; r) = Id, 


where A* is defined by (2.4), which we reproduce here: 

-4*(t) = {Q{A- fjL Id)(5“^) xo + x,C),C*(e^^^T x,C))- 

Recall that h = 2j3 and C = 1/3 here. The change of basis Q is given by Lemma 5.3 


the real part of the branching eigenvalues ^ is introduced in Lemma 5.1, and (x*, C*) is the 
bicharacteristic flow, solving (2.5), which we reproduce here: 

(5.7) 9fX* = Xo + e^/^x*, C*), 5tC* = xq + C*)- 

The block decomposition of A given by Lemma |5.3| induces a block decomposition of A*. 
We focus on the top left block A(o)* in A*. As per Lemma 5.3, its bottom left entry involves 
the function 

(5.8) r**(e,t,x,C) := r*(xo + e^/^x*(e^/^t,x,C),C*(e^'^^La^,C) )■ 

The bicharacteristic flow ( |5.7| ) satisfies 

(5.9) Xo + e^/^x*(e^/^L x, 0 = 2:0 + + 0{et), x, ^) = .^ + 0{e), 

uniformly in (x, with |x| + |^—^o| < Here notation 0(e) refers to a uniform bound of form 
< . In particular for t bounded from above by a power of | lne|, there holds 0(et) = 0(e). 
Thus r** defined in (5.8) satisfies 

T** = 6»*(e^/\,,f) + 0 (e), Oi,(x,C) := r*(xo + x,^). 


By (5.4), we see that 0* as defined above satisfies the conditions of equation (2.1). In 


accordance with ( |2.1[ ), we let t*(e,t, x,^) := e“^/^0*(e^/^x, ^). There holds 

(5.10) - r**) = e^^^(t - ^)) + 0(e^/^) = e^^^(t - U) + 0(e^/^), 

uniformly in (x,.^) with |x| + |C ~ Col < 
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5.3 


□ 


Lemma 5.4. The top left block ^(o)* ofA^ in the block decomposition of Lemma 5.3 satisfies 

(“)* V - *.)/o + 0 ) ' 

with notation /o(e,x,^) = eo(0, xq + e^/^x,^), so that /o > 0 for {x,f^) near (0,^o)- 

Proof. The bottom left entry of involves the function eg = e(/i) evaluated in the 

time-rescaled and advected frame. In view of ( |5.9[ ), there holds by regularity of e and p, 

eo(e^''^t,xo + = eo(0,xo + 0(e^/^). 

The bottom left entry of also involves t — Ti, and (t — in the time-rescaled 

and advected frame. In view of (5.10), these functions contribute to 

e“i/3(e^/3t - r**) = e^/^{t - f*) -h 0{e^/^) 

and — = 0(e^/^). We may conclude with Lemma 

By Lemma 5.4 the flow 5(o) of solves 

(5.11) + _,i/3( 

where C := C(e,f,x,^ ^ 0 ) ’ ^ 1 for 0 < f < T{e) and |x|-|-|^—^o| < d. 

The coefficient /o = /(e, x, ff) satisfies /o > 0 for |x| -|- |^ — Col < <5- 

5.4. Reduction to a perturbed Airy equation. Let 

Zl(x,0 J). 

SO that D is well-defined and invertible on |x| < d, |C — Col < and 

(5.12) Z{T]t) ■.= DS{o){fQ^^^T + U-,fQ^^^t + u), 

where /o, f*, D, S'(o) and Z all depend on (x,C)- For future use, we note that 

-*(e/o) ^^32(12 


.S(o)=e2/3^S(o), 5(o)(T;r) = Id, 


D{x,0 


Zu Zi2 
Z21 Z22 


zi(x,c)-' = 


^11 


*(e/o) 


- 1/3 


Z21 


Z22 


(5.13) 
and 

(5.14) Il(x,C)-3 
Lemma 5.5. On |x| < S, |C — Col < d, the map Z satisfies the perturbed Airy equation 


) 0(31.0 = ( Ai/3 

^21 Z22 J ^ V -ii£for^Z2l Z22 


0 1 
t 0 


Z = £^I^CZ, Z{t]t) = D, 


(5.15) Z'+ 

where C ;= {DCD~^){fQ ^^^t + t.^), with C as in (5.11) 
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Proof. It suffices to use (5.13) and observe that, by Lemma 5.4 


□ 


From the above, we will deduce lower and upper bounds for *S'(g), by comparison with the 
vector Airy function Z, defined as the solution of 


(5.16) 


Z' + 


0 1 
t 0 


Z = 0, 


Z(r; r) = Id. 


5.5. Bounds for the Airy function. We will use (5.15) to show that the symbolic flow 
grows in time like the Airy function, for which the following is known (see for instance [5], 
chapter 7.6). 

Lemma 5.6 (Airy equation). Let Ai he the inverse Fourier transform of , and j = 
^i-kIz^ Tfie Junctions Ai, Ai(j-) form a basis of solutions of the ordinary differential equation 
y” = fy, and there holds 

Ai(t) = * + 0(t-s/2)), t ^ +«,, 

Iff'K 


Ai(it) = 


20F 


^/6e(2/3)*^/^t-l/4(l + 0(t-3/2)), ^ ^ 


Ai(-jt) = + 0(t-3/2)), t ^ +oo. 

2y TT 

From the above Lemma, we deduce uniform bounds for the time derivative Ai' : 
g(2/3)i3/2|^-/(^)| ^g-(2/3)t3/2|A-/(^-^)| ^ |Ai'(-t)| + |Ai'(-jf)| < C{l+tf/\ 

for some C > 0, for all f > 0. By Lemma 5.6, the solution to (5.16|) is 


Z(r;t) = 


-jAi'(j'r)Ai(t) + Ai'(r)Ai(jt) -Ai(jr)Ai(t) + Ai(r)Ai(jf) 


W{t) V 4 Ai'(jr)Ai'(f) - j'Ai'(r)Ai'(jf) Ai(jr)Ai'(t) - j'Ai(r)Ai'(jt) 
where VF is the Wronskian, satisfying 

IF(r) := Ai(jr)Ai'(r) - j'Ai'(jr)Ai(r) = ^(-\/3 + z). 

The bounds for Ai and Ai' imply the upper bound, for 0 < r < f : 

(5.17) |Z(r;f)| < (7(1 + |r|)3/'^(l + |f|)^/^eAi(T; f), 

and the lower bound 


(5.18) 


( 1 0 )Z(0;t) 


> ceAi(0;t), 
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for some c > 0 independent of T,t, where the growth function eAi is defined by 

(5.19) eAi(r;t) = exp := max(x, 0). 

We note that eAi is multiplicative: 

(5.20) eAi(r; t')eAi{t'] t) = eAiir; t), for all r, t', t'. 

Remark 5.7. If we had assumed d\PdtP < 0, then we would have had to eonsider the Airy 
eondition for negative times. Lemma ^ would then have yielded polynomial bounds for the 
symbolie flow. 

5.6. Bounds for the symbolic flowf. Let 

(5.21) 0(L X, 0 := fo{x, - tfle, x, f,)). 

Our goal is to verify the bounds of Assumption 


2.1 


for S'(o) in the elliptic domain V defined 


in (2.2). We reproduce here the definition of P : 

P := {(r;t,x,0, U{e,x,() < t < t < T{e), |x| < 5, |^ - ?o| < 


Lemma 5.8. There holds the bounds in domain P : 

1 e-V3 


D-^Z{e{T);Q{t))D\ 


M/3 


1 


eAi(0('r);0(t)), 


with CAi defined in (5.19). 


Above, < means entry-wise inequality “modulo constants”, as defined in (2.12). 


Proof. There holds 0 < 0(r) < 0(t) < 0(T(e)) in domain P. Bound (5.17) states that 
there holds |Z(0)| < eAi(0). Then (5.14) implies the result. □ 


From Lemma 5.8 we now derive bounds for -5(0). Given that 5(o) is expressed in terms 
of Z (5.12) and that Z is a perturbation of Z, we find ourselves in a situation very much 
like the one encountered in Section 3.12[ Accordingly, the proof of the following Corollary 
borrows from Section [3.12[ in particular from the proof of Corollary 3.9 and Corollary |3.8[ 


Corollary 5.9. The flow 5(o) of the top left block A*(o) in A*, solution of (5.11), satisfies 
the bounds: 


(5.22) 


|5(o)l 


1 £-V3 

eV3 1 


eAi(0), 


Proof. By Lemma 5.5 and definition of Z (5.16), there holds 

Z(r;t) = Z(r;t)L> + e^/3 [ Z{t']t)C{t')Z{Tfl') dt'. 


By definition of Z (5.12), there holds 5(o) = D ^Z(0). Thus 

pe{t) 

Sfo){T;t) = D~^Z{e{T);e{t))D + / D-^Z{t';e{t))C{t')Z{e{Ty,t') dt'. 

Je{T) 
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5.5 


to be equal to [D ^CZ))(/q +1*), we obtain 


Since C is defined in Lemma 

S^0){r-,t) = D-^Z{e{T)-e{t))D 

Je{T) 

_1 /o 

The change of variable t' = O(t'), corresponding to t' = + transforms the above 

integral into 

^2/3^i/ 3 0(t))Z)C(r')5(o)(r; r') dr'. 

We now factor out the expected growth in view of applying Gronwall’s lemma, as we did 
before in the proof of Corollary 3.9 we let 

Sfo) :=eAi(0)-'5(o), and Z^(r'; t) := eAi(0)-^d?"^Z(0(r'); 0(t))i?. 

By the multiplicative property ( 5.20[ ) of the growth function cai, we find 

Slo)iT-,t) = z\T;t) + e‘^^^fo f Z^(r'; t)C'(r')5(o)(r'; t) dr'. 


We now rescale the top right and bottom left entries, as we consider the equation in 5 


fo)’ 


with notation introduced just above (3.56) in the proof of Corollary 3.8 In view of (3.56), 
there holds 

(5.23) SUr-,t) = Z\r-t)+e^/^fo f Z\r';t){e^/^C{r'))SUT';t) dr'. 


5.8 


There holds e^/^C(t) = 0(t), and t is bounded by some power of | lne| in T>. Lemma 
implies that |Z^| < 1. Hence Gronwall’s lemma implies the bound liS/gHr;t)| < 1, which 


corresponds precisely to (5.22). 


□ 


Lemma 5.10. There holds the lower bound 

0 
1 


( 1 0)S(0)( 

for some universal constant cq > 0. 


>coe 3/3g^.('0)^ 


Proof. Consider representation (5.23). We focus on the top right entry. The lower bound 
(5.18) for the vector Airy function states that the top right entry of Z is bounded from 
below by eAi. By (5.14), this implies 

(5.24) (10) D-^ZD 


> Coe 3/3eAi, 


for some cq > 0 independent of r, t. Borrowing notation from the proof of Corollary |5.9[ 
this means that the top righ t entry of Z is bounded away from zero, uniformly in time. 

that |*S(^o)l ~ 1 ^^0 \Zi \ < 1. Thus from (5.23) and (5.24) we 
in V. □ 


5.9 


We know from Corollary 
deduce the result, since t^ 
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We observe that there holds, for CAi defined in (|5.19) and 0 defined in (5.21): 


eAi(0) = e, 


• 7 ? 


with 7 (x,^) ■■= U = £ + 


Lemma 


5.10 


being equal to the constant vector 


2.1 


for 5, 


where r* is given the implicit function theorem in Section |5.1 
verify the bounds of Assumption 

0 
1 


( 0 ), 


Hence Corollary 5.9 and 
with 7 ’'“ = 7 “ = 7 , and with e 


In order to complete the verification of Assumption 2.1 and thus conclude the proof of 


Theorem 1.3, it only remains to show that the other components of the symbolic flow do not 


grow faster than ^o)- Tbis follows directly from the simplicity hypothesis in Theorem 1.3 


Indeed, by the simplicity hypothesis, we may smoothly diagonalize the other component 
of A* near (0, xo,^o) (use for instance Corollary 2.2 in [25]). The eigenvalues of 
are real near (0,xo,Co)- The equation for the symbolic flow of splits into scalar 

differential equations, with purely imaginary coefficients. Thus the symbolic flow of 
is bounded. 


6. Proof of Theorem II. 6t smooth defect of hyperbolicity 


It suffices to verify that, under the assumptions of Theorem |1.6[ Assumption 2.1 holds 
with parameters 


£ = 1 , h = l/2, C = 0 , ;U = 3f?eA±, f* = 0 , 


where A± are the bifurcating eigenvalues, as given by Proposition 1.4 


6.1. Block decomposition. As in the proof of Theorem 1.2 we may smoothly block 
diagonalize A by a change of basis for small t and {x,S,) close to (xo,Co)- Then 

identity (4.2) holds, and we focus on block A(o), of size two, such that 

spA(o)(0,xo,Co) = {Ao}, Ao G M, 


where (xo,^ 0 )Ao) are the coordinates of wq G P which intervenes in Hypothesis 1.5 
Hypothesis and Proposition 


1.4 


By 


the eigenvalues A± of A(o) branch out of the real axis 
at f = 0, for all (x,^) in a neighborhood of {xo,^q). We define fj. to be the real part of these 
eigenvalues. The corresponding equation for the symbolic flow is 

(6.1) dt5(o)+e"^/^(A(o) - |u)(e^/7,xo+ =0, 

where (x*,^*) are the bicharacteristics of fi. 


6.2. Time regularity and cancellation. By Proposition 1.4, the eigenvalues X± are dif¬ 
ferentiable in time, at f = 0 and for all (x,^) near (xo,Co)- Indeed, Hypothesis 1.5 implies 
that conditions ( 1 . 12 ) are satisfied in a whole neighborhood of (xo,^o)- We may thus write 

(6.2) A±(e^/^f,xo + e^/^x*,^*) - /i(0,xo + e^/^x*,^*) = i£^^‘^tX±{e,t,x,^) + o(e^/^), 

uniformly in f = 0{\ lne|) and (x,^) near (xo,^o )5 where (x*,^*) is evaluated at (e^/^t,x,^), 
and where 


(6.3) 


A±(e,0,x,O = A±(0,xo + e^/^x*(0,x,0,^*(0,x,0) G 
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Consider the 2x2 matrix It has one semi-simple eigenvalue /r(0, x,^) (the 

assumption of semisimplicity is part of Hypothesis (1.5)). Thus 

^(O)(0,a:,0 = 

In particular, by regularity of the entries of A, 

(6.4) + = ti(o)(e, 0, x, 0 + t, x, 0, 

where B is uniformly bounded for e close to 0, t = 0(|lne|*) and (x,^) close to (xqj'^o)- 
Thus equation (6.1) takes the form 

(6.5) dtS(o) + tA(o) (e, 0, x, OS'(o) = t, x, 05'(o) ■ 


The key cancellation that takes place in (6.4) transformed equation in 5(0) into an au¬ 
tonomous equation with a small, linear, time-dependent perturbation. The eigenval ues of 


ho) 


are A±(e, 0,x,.^) from (6.2)-(6.3). These eigenvalues are distinct by Proposition 


1.4 


6.3. Bounds for the symbolic flow. The solution S to 

dtS -I- itH(o)(e, 0,x,^)5 = 0, 5(r;r) = Id 

is 

5(r; t) = exp (^-ii(o)(e, 0, x, 0(i^ - t^)/2 ) • 

The eigenvalues of H(o), being distinct, are smooth in (e,x,^) (see for instance Corollary 
2.2 in [25|1. In particular, there holds 

A±(e,0,x,O = A±(0,0,x,O + 0{e) = 9mdtA±(0,xo,C*(0,x,^)) 

locally uniformly in Let A+ be the eigenvalue with positive imaginary part, and 


7(a;,0 := ^A±(0,0,x,0 = (9tA+(0 ,xq,^*( 0,x,0)- 

Then, 

(6.6) |5(r;t,x,0| < exp - r^)), 

and, since H(o) is smoothly diagonalizable, for some smoothly varying vector e(x, .^) there 
holds 


(6.7) 


|5(T;t,x,^)e(x,0 I > exp {'y{x,^){t‘^ - r^)). 


Perturbation arguments already encountered in Section |3.12 (specifically, in the proof of 
Corollary |3.9| ) show that the bounds (6.6)-(6.7) for S yield similar bounds for the symbolic 
flow 5(0) solution to (6.5). These bounds verify the upper and lower bound (2.11) and (2.10) 
from Assumption |2 . l] 


For the other components of the flow, we use the simplicity assumption in Theorem 1.6 


as we did in the last paragraph of Section 5.6 in the proof of Theorem 1.3 
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7. Examples 

7.1. One-dimensional Burgers systems. The 2x2, one-dimensional Bnrgers system 


(7.1) 




where F and h are smooth and real-valued, has a complex structure if h is constant. In 
the case 6 = 1, F = (0,1), a strong instability result for the Cauchy-Kovalevskaya solution 
issued from (n?, 0 ), where is analytic and real-valued, was proved in [ 12 ] • 

We assume 6 > 0, and the existence of a local smooth solution (p = (^i, (p 2 )- The principal 
symbol is 

4>2 4>i 


A{t,x,P,) = ^ 

Without loss of generality, we let ^ = 1. The eigenvalues and eigenvectors are 

1 


X±= 4>i± i4>2bi^), 
The characteristic polynomial is 


e± = 


(1 + 6(,/.)2)i/2 


±ib{(l)) 

1 


P= {X - pif + bipfpl 

Initial ellipticity. If 4>2iO,xo) / 0 for some xq G M, then the principal symbol is elliptic 
at t = 0, and Theorem appplies. 

Smooth defect of hyperbolicity. Consider the case ^ 2 ( 0 , x) = 0. We cannot observe a 
defect of hyperbolicity as in Theorem 1.3 since the eigenvalues are smooth in time. Via 
Proposition |1.4[ we see that Theorem 1 1. 6 | holds as soon as 


(7.2) F 2 ((?!)( 0 , xo)) / 0, for some xq G M. 

In the case h{u) = b{u 2 ), then is a system of conservation laws 
dtui + dxfiiu) = Fi{u), dtU 2 + dxf 2 {u) = F 2 {u), 

with fluxes 

2 f'^2 

h{u) = -uf - yb{yf dy, f2{u) = U1U2. 

If, for instance, F{u) = (0, u\) and 6 (^ 2 ) = 1 -|- then the system is ill-posed for all data. 
7.2. Two-dimensional Burgers systems. Consider the family of 2 x 2 systems in M2 ; 

uidx, -b{uYu2{dx2+dxf) 


dtu + 


U2{dxi T dx2^ 


u = F{u). 


■'X\ ^ '-'X2) ^idxi 

We assume 6 > 0, and the existence of a local smooth solution (p = (pi, p 2 )- The principal 
symbol is 

601 -(6+6)&(0)^02 
(6 + 6)02 601 


A = 
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The eigenvalues and eigenvectors are 

=t *(6 + C2)4>2b{4>), 


e± 


1 

(1 + fe (</.) 2 ) l /2 


±ib{(j)) ^ 


Initial ellipticity. If (/>2(0,xo) / 0 for some xq G 
elliptic at any (^ 1 ,^ 2 ) £ such that ii+ ^ 2 ^ 0. 


then the principal symbol is initially 


Smooth defect of hyperbolicity. Consider the case (/)2(0, x) = 0. By Proposition 1.4 the 
assumptions of Theorem 1.6 are satisfied under condition (7.2). 


7.3. Van der Waals gas dynamics. The compressible Euler equations in one space di¬ 
mension, in Lagrangian coordinates are 


J dtui + dxU2 = 0, 

\ dtU2 + dxp{ui) = 0 . 

We assume that the smooth pressure law p satisfies the Van der Waals condition 


p'(ui) < 0, for some ui G M, 

and assume existence of a smooth solution </> = {(j)i,4>2)- The principal symbol at ^ = 1 is 


The eigenvalues are 


A = 


0 1 
P'(f>i) 0 


A± = (p'iM)^''^- 


Initial ellipticity. If p'((^i(0, xq)) < 0 for some xq £ M, then Theorem 1 1. 2| applies. 

Non-semi-simple defect of hyperbolicity. If p'((/)i(0, x)) > 0 for all x (initial hyperbolicity) 
and p'{(fi{0,xo)) = 0 for some xq (coalescence of two eigenvalues), if 


p"XQ))dx(p2{0, xq) > 0, 


then condition (1.11) holds, and Theorem |1.3| applies. 


7.4. Klein-Gordon-Zakharov systems. Consider the family of systems in one space 
dimension 


(7.3) 



V 

u 

m 

n 






indexed by a G M, c G M\{ —1,1}. We assume existence of a smooth solution cj) = (u, x, n, m). 
The principal symbol at ^ = 1 is 


/ 

0 

1 

a 

0 \ 


1 

0 

0 

0 


a 

0 

0 

c 

V 

—2u 

-2v 

c 

0 / 


(7.4) 
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The case a = 0. The principal symbol is block diagonal, and there are four distinct 


eigenvalues {±1, ±c}. This implies that (7.3) is strictly hyperbolic, hence locally well-posed 
in H^, for s > 3/2 (see for instance Theorem 7.3.3, [E]). It was observed in [ 2 ] that for 
c ^ {—1,1} and a = 0, system (7.3) is conjugated to a semi-linear system, which implies a 
sharper existence result: 


Proposition 7.1 ([2|, Section 2.2). If c ^ |—1,1| and a = 0, the system (7.3) is locally 
well-posed in for s > 1 / 2 . 


Proof. The change of variables 

{u, v) = {u + v,u — v), 


(h, rh) = ( n -|- m — 


1 


1 — c 




-V , n — m — 




1 — c 


transforms (7.3) into the system in U := {u,v,h,rh) : 

/ 0 1 0 0 \ / 


(7.5) 


dtU + 


10 0 0 

0 0 0 c 

y 0 0 c 0 y 


dxU = (n -M) 


— V 

u 

— 2(1 — c)^"uv 
\ - 2 (l + c)-i 




uv 




for s > 1/2. 


System (7.5), being symmetric hyperbolic and semi-linear, is locally well-posed in H‘ 


□ 


The case a 0. By Proposition 7.1, system (7.3) takes the form of a symmetric pertur¬ 


bation of a well-posed system. The characteristic polynomial of the principal symbol (7.4) 
at ^ = 1 is 


P{t, X, A) = (A^ — c^)(A^ “ 1) “ o^A^ -|- 2q;c(u -|- uX). 

Consider an initial datum for (u(0),u(0),n(0),m(0)) such that, for some xq £ M, 

(f 

(7.6) tt(0,xo) = 0, v{0,xo) = - — , acdxu{0,xo) > 0. 

2a 

The first two conditions in ( |7.6[ ) imply that at wq = (iCO) 1)0) there holds 

P(0,wo) = dxPiO,uJo) = 0. 

The third condition in ( |7.6[ ) implies 

{dtPdlP){0,u}o) = {2acdtv{0,xo)){-l - (? - o?) = (2acdxu{0,xo)){l -\-c^ -\- a^) > 0, 

so that the third condition in (7.6) implies condition (mj- Theorem |1.3| thus asserts 
instability of the Cauchy problem for (7.3) in the vicinity of any smooth solution cj) satisfying 
(7.6) at t = 0. 

In particular, for any given ao > 0, we can hnd initial data, depending on uq, such that 


with a = 0 is well-posed whereas (7.3) with a = ckq is ill-posed. Such initial data are 
O(l/ao) in L° 
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Appendix A. Proof of Proposition 11.41 


The principal symbol can be block diagonalized, with a 2 x 2 block Aq with double real 
eigenvalue Aq at ( 0 ,x,^), and an {N — 2 ) x {N — 2 ) block which does not admit Aq as an 
eigenvalue at (0, x,^). Throughout this proof {x,^) are fixed and omitted in the arguments. 
The characteristic polynomial of A factorizes into P = PqPi, where Pi(0,a;o) 7 ^ 0 and Pq, Pi 
have real coefficients. We may concentrate on Pq : 

Po{X) = A^ — Atr Aq + det Aq. 


The eigenvalues X± of Aq at are 

(A.l) A±(t) = ^tr Ao(f) ± A(f) := (tr Aq)^ - 4det Aq. 


By assumption, these eigenvalues coalesce at t = 0, so that A(0) = 0. The goal is then to 
prove equivalence ( 1 . 12 ). 


If the left proposition in (1.12) holds, then A{t) = —aP + 0{P), with a > 0. Thus 
dtA{0) = 0 ; on the other hand. 


clfA(O) = 2tr Ao(0)clttr Ao(0) — 49* det Ao(0) = 4Ao9ttr Ao(0) — 49tdet Ao(0) 
= -4(9*Po)(0). 


Besides, 9|A(0) < 0 ; on the other hand, 

dfA(O) = 2(9*tr Ao(0))^ + 2tr Ao(0)9^^tr Ao(0) — 49^^ det Ao(0), 

implying, since tr Ao(0) = 2Ao, 

9*^A(0) = 2(9*tr Ao(0))^ + 4Ao9^tr Ao(0) — 49*^ det Ao(0) 

= 2(9*9aPo)' - 2d!Pod^Po(0), 

which gives indeed (9^\P)^ < 9|P9|P at t = 0 . 

The converse implication is proved in the same way: the right proposition in (1.12) 
implies 9*A(0) = 0 , 9f A( 0 ) < 0 , as shown above, and this implies that the eigenvalues in 
(A.l) are differentiable and leave the real axis at t = 0. 


Appendix B. Symbols and operators 
Pseudo-differential operators in e^-semi-classical quantization are defined by 

(B.l) opi;{a)u := f P^'^a{x,e^^)u{^) 0 < e, 0 < h. 

Here h = 1 /( 1 + £), as in Assumption [ 2 ^ Above a is a classical symbol of order m: a ^ 5™, 
for some m G M, that is a smooth map in (x, ^), with values in a finite-dimensional space, 
such that 

(B.2) ||a||^_^ := sup 01 < 00 , 

|a|<r,|/3|<r 
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The family || • of e-dependent norms is defined by 

Wk. »£R. 0:=(1 + I'l")'k 

Introducing dilations {d^) such that {deu){x) = we observe that there holds 

(B.3) ||dert||H® = op£(a) = (i“^op(a)d£, d{x, ^) := a{e^x, . 

Proposition B.l. Given m G M, a G 5"*, there holds the bound 
(B.4) \\op^{a)u\\L 2 < \\a\\^^c{d)\Me-m, 

for all u G H~'^, for some C{d) > 0 depending only on d. If m = 0, there holds the bound 

(B.5) l|oPe(a)u||i2 < ^ supJ5X-,0lLi(Rd)||^^||L2. 

0<|a|<d-|-l^ 


Proof. By use of dilations (B.3), we observe that op£(o)u = opi{{f)~^d){D)^d^u. Bound 
(B.4) with any C{d) > [d/2] -|- 1 then follows for instance from Theorem 1.1.4 and its proof 
from m- Bound ( |B.5[ ) is proved in Theorem 18.8.1 from volume 3 of [5]. □ 

Proposition B.2. Given ai G 02 G re G N, 

op£(ai)op£(a2) = ^ e'*‘'op£(aitlga2)-he'*(”+^^op£(i?n-ri(ai,a2)), 

0<q<n 


where 

(B.6) 


aittga 2 = ^ d^aid^a2, 


W\=q 


and i2„+i(oi,02) G (*^+1) satisfies 


|op^(i?n+i(ai,a 2 ))re|| 2,2 


< 


; llmi,C(d) 11^3, 0.211 m 2 ,C(d) Il'^l|£,mi-|-m2—n—1 j 


with C{d) > 0 depending only on d, for all u G n i_ 

Proof. Based for instance on Theorem 1.1.20, Lemma 4.1.2 and Remark 4.1.4 of m, and 
the use of dilations (B.3). □ 

Specializing to symbols with a slow x-dependence, we obtain: 

Proposition B.3. Given ai G 02 G , if 0.2 depends on x through e^~^x, there 

holds 

||(op£(ai)op£(a2) - opg(aia2))re||£,s < e||ai|| 

mi,C{d) 11^2||'m, 2 ,C(cZ) ||^||£:,s+mi+m 2 —1? 
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Appendix C. On extending locally defined symbols 


Our assumptions are local in (x, .^) around {xq, ^q). Accordingly the symbols Q and jjL that 
appear in Assumption 2.1 are defined (after a change of spatial frame) only locally around 
(0,^o)- We explain here how to extend the locally defined family of invertible matrices 
Q{x, into an element of with an inverse (in the sense of matrices) which belongs to 

50 . 


The spectrum of <5(0, ^ 0 ) is a finite subset of C. In particular, we can find a G M such 
that the spectrum of <5(0, .^ 0 ) is included in C\e*“M_. By continuity of the spectrum, for all 
(x, ^) close enough to (0, ^ 0 )) the spectrum of Q(x,^) does not intersect the half-line 
Let 5 > 0 such that this property holds true over Bs = B{0,S) x B{^q,6). We may then 
define the logarithm of matrix iii Bs by 


Log(e-*“g)= [ (e-*"Q-Id)((l-t)Id+te-*“g) ^ dt, 
Jo 

and the notation Log is justified by the identity 


(C.l) expLog(e *"g) = e *“g, in Bs- 

Let <j(x,^) be a smooth cut-off in C'“(M'^ x M*^), such that 0 < cr{x,^) < 1, with u = 1 on 
a neighborhood of (xo,^o)) and such that the support of a is included in Bs/ 2 - Let 

R{x, 0 = cr{x, 0Log(e“*"g(x, 0) + (1 - o-{x, 0) W, in Bs- 

We may extend smoothly ii by i? = Id on the complement of Bs in Then for all 
(x, ^) G the matrix 

g(x,0 = expi?(x,0 

is smooth and invertible. There holds 


inf det g > 0. 

R 2 d 


Indeed, the infimum over the closed ball Bs is positive, by compactness and continuity, and 
the determinant is constant equal to outside Bs- Thus the norms |g(x,^)| and |g(x,^)“^| 
are globally bounded over Since Q is constant outside a compact, this implies Q £ S^, 
and Q~^ G S^. Finally, by (C.l) and definition of the cut-off a, there holds 


g(x,0 = e *“g(x,|), for (x,,^) close to (0, ^o)- 
Thus e*“g is an appropriate extension of Q. 


Appendix D. An integral representation formula 


We adapt to the present context an integral representation formula introduced in |24j . 
Consider the initial value problem, posed in time interval [0,r(e)], with the limiting time 

T(e) := (T*| loge|)^'^^^'''^\ for some T* > 0, some i >0 : 

(D.l) dtu + op^{A)u = g, u(0) = uq, 


where A = A{e,t) belongs to for all e > 0 and all t < T{e). Recall that opg(-) denotes 
e^-semiclassical quantization of operators, as defined in (B.l). The parameter h belongs 
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to (0,1]. The datum uq belongs to L^, and the source g is given in T(e)],L^(]R'^)). 

Denote Sq the flow of —A, defined for 0 < r < t < T(e) by 


dtSoir-1) + ASo{t-, t) = 0, 5o(t; r) = Id. 

For some G N*, denote {*S'q}i<q<go the solution to the triangular system of linear ordinary 
differential equations 


(D.2) 


dtSq + ASq + ^ A%^Sq^ = 0 , Sq{T\ t) = 0 , 


<?i+g 2 =g 

0<gi 


with notation Ug introduced in (B. 6 ). 

Assumption D.l. The symbol A is compactly supported in x, uniformly in e, t, and there 
holds bounds 

and 

\d^d^^Sq{e, r; t, xA)\ < exp . 

for some 0 < ( < h, some 7 > 0, for all x, all t <T{e). 

Denote S := s^^Sq. Then opj(S) is an approximate solution operator for (D.l): 

o<g<go 

Lemma D.2. Under Assumption D.l[ if go is large enough, depending on (, h, 7 and T*, 
there holds the bound 


(D.3) ^toPeCS) + op^(A)op^(E) = p, 

where p satisfies for 0 < t <t <T{e), for all u G 
(D.4) l|p(^;i)^^llL 2 < e||w||L 2 - 

Proof. By Proposition |B.2[ 

op^(A)op^(Sq) = op^(A5g) + e'^'^op^{A]\q>Sq) + e(®+^)^op^(iigo+i(A, Sq)), 

^<g'<go 

and, summing over 0 < q < qo '■ 

op,(A)op,(E) = op,(AS) + 

0<g2<go 

^<gi<go 


where R := 


e'^^op^{Rqq+i{A, Sq)). Besides, by definition of the correctors (D.2), 


o<g<qo 

-(9top£(S) = op^(AS) + Y 


^<gi+g2<qo 

0<qi 
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Comparing with the above, we find that (|D.3[) holds, with 


-P ■= 




9o+l<qi+92<2go 

0 <q 2 <qo 

Since A is compactly supported in x, so are the correctors S'^, for g > 1, and the derivatives 
of So- Thus under Assumption |D . 1 there holds 

\dx{MgiSg 2 )\Li(Kd) ^ exp 

uniformly in Hence, by Proposition |B.1[ the bound 

||0p,(Attgi5g2)|L2^i2(Krf) < OXp - 

Thus a control of the —)• norm of the first term in p by 

^{qi+q2)h^-Cil+qi+q2) g^p (jAY < £(^“C(90 + 1)-C-7A^ 

qo+i<qi+q 2 < 2 qo 

0<'?2<IJ0 

over the interval [0,T(e)], implying the desired bound as soon as 

{h - C){qo + 1) > 1 + C + jT^- 

Besides, by Proposition |B.2[ there holds 

||oPe(-Rqo+l(-^>'S'g))|| £, 2_^^2 ^ ||9| •^llo,C(d) 'S'(j||o,C(d)) 
and by Assumption |D.1[ 

l|9f+U||c(d) < e-(®+i+cW)C, 115^0+15^11^^^^ < ,-ai+q+cid)) g^p 

This gives a control of the —)• norm of the second term in p by 

g-lw+l)^ y ^ ^qh^—(qo+l+C{d))(^—{l+q+C(d))Q g^p 
0<g<<Jo 


We conclude that (D.4) holds if qo satisfies 

— C)(Q'o + 1) > 1 + 7 T* + 2C{d)(, 


which can be achieved since in Assumption D.l we postulated <h. 


□ 


Theorem D.3. Under Assumption D.l, the initial value problem (D.l) has a unique solu¬ 
tion u G C'°([0,T(e)],L^(]R'^)), given by 


(D.5) 


u = op£(S(0; t))no + / op£(S(t'; t))(Id+£i?i(t'))( 5 '(t') + eR 2 {t')uo) dt', 


where Ri and R 2 are bounded: for all u G L^, 

(D.6) \\Ri{t)v\\L2 + \\R2{t)v\\L2 < \\v\\l2, 

uniformly in e and t G [ 0 , T{e)]. 
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Proof. Let h G -L°°([0, T(e)], By Lemma D.2, the map u defined by 

It := op£(S(0; t))ito + [ dt' 


solves (D.l) if and only if, for all t, there holds 

(D.7) {{ld+po)h){t) = g - p{0;t)uo, 

where po is the linear integral operator 


Po 


f G C'°([0, r(e)], L^) —>■ —)• f p{T]t)v{T) dr^ G C^{[0,T{e)], L"^). 


By (D.4), there holds 


sup ||(po^')(^)||I ,2 < e sup \\v{t)\\L 2 . 

0<t<T{£) 0<t<r(e) 


Thus Id+/?o is invertible in the Banach algebra of linear bounded operators acting on 
C'*^([0, r(e)], This provides a solution h to (D.7), and we obtain representation 

(D.5) with Ri := e“^((IdTpo)”^ ~ Id) and R 2 = —/o(0; •). Bound (D.6) is a consequence of 
(D.4). Uniqueness follows from Cauchy-Lipschitz, since ^ G 5®. □ 
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